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CHAPTER  1 
INTRODUCTION 

The  objective  of  this  thesis  is  to  analyze  the  performance  of  asynchronous  binary  PSK  direct- 
sequence  spread-spectrum  multiple- access  (DS/SSMA)  communications  over  an  additive  impulsive 
noise  channel.  In  recent  years,  there  have  been  a  number  of  studies  of  such  communication  systems 
with  the  additive  white  Gaussian  assumption  for  channel  noise.  In  particular,  many  researchers  have 
devoted  their  efforts  to  finding  efficient  methods  of  obtaining  approximations  and  bounds  for  the 
average  bit-error-probability  of  the  conventional  correlation  receiver  within  this  model  [11-12],  [26- 
27],  and  [43-44].  Pursley  [26]  developed  an  approximation  for  the  average  error  probability  based  on 
the  average  signal-to-noise  ratio.  Using  moment-space  ideas,  Yao  [44]  developed  upper  and  lower 
bounds  on  the  error  probability.  Later,  Wu  and  Neuhoff  [43]  used  series  expansions  and  Gauss  qua¬ 
drature  rules  to  obtain  estimates  for  the  performance  of  the  conventional  linear  correlation  receiver. 
Recent  contributions  to  this  problem  include  the  methods  of  approximations  given  by  Geraniotis  and 
Pursley  in  [1 1]  and  [12],  which  are  based  on  series  expansions  and  the  integration  of  the  characteristic 
function,  respectively.  Each  of  these  methods  has  its  advantages  and  disadvantages,  and  the  choice  of 
method  for  a  given  application  ultimately  depends  on  the  system  parameters  and  required  accuracy. 
Taken  collectively,  a  number  of  efficient  techniques  are  available  for  evaluating  the  performance  of 
the  linear  correlation  receiver  in  additive  white  Gaussian  noise  in  the  DS/SSMA  environment 
Among  the  many  other  contributions  to  the  performance  analysis  of  DS/SSMA  communications  are 
the  studies  of  Anderson  and  Wintz  and  others  (  see  [5] )  who  have  analyzed  DS/SSMA  systems  with 
a  hard  limiter  in  the  structure  of  the  receiver,  and  of  Verdu  (  [41]  and  [42]  )  who  has  considered  an 
optimum  detection  algorithm  that  yields  a  minimum  error  probability  DS/SSMA  receiver. 

One  common  assumption  among  the  contributions  referenced  above  is  that  of  additive  white 
Gaussian  noise.  Although  this  assumption  is  quite  appropriate  for  many  applications,  it  is  well- 
known  that  many  noise  environments  arising  in  practice  are  poorly  modeled  by  Gaussian  statistics. 
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In  particular,  man-made  electromagnetic  interference  ( or  noise  ),  and  a  great  deal  of  natural  interfer¬ 
ence  as  well,  is  basically  "impulsive"  in  nature  ;  i.e. ,  it  has  a  highly  structured  form,  characterized  by 
significant  probabilities  of  large  interference  levels.  This  is  in  contrast  to  the  more  "entropic"  Gaus¬ 
sian  noise  processes  inherent  in  transmitting  and  receiving  elements.  This  impulsive  or  structured 
character  of  the  interference  can  significantly  degrade  the  performance  of  conventional  (  linear  ) 
demodulation  systems,  which  are  designed  to  operate  most  effectively  against  the  commonly 
assumed  Gaussian  background  noise  processes.  However,  by  proper  system  design,  this  noise  struc¬ 
ture  can  often  be  exploited  to  yield  better  performance  than  would  be  obtained  in  the  Gaussian  case. 
There  have  been  various  efforts  over  the  past  two  decades  in  the  area  of  signal  detection  in  impulsive 
noise,  primarily  within  the  context  of  single-user  channels  [3],  [7-10],  [13],  [18-23],  [36-38],  and 
[40],  Most  of  these  studies  are  concerned  with  the  problem  of  detecting  the  presence  of  a  weak  signal 
in  a  noisy  environment  [3],  [7],  [14],  [18],  [21-23],  [36-38],  and  [40],  and  they  illustrate  that  models 
for  the  first-and  second-order  probability  distributions  of  the  noise  are  usually  necessary  in  order  to 
achieve  good  performance.  Since  optimum  systems  for  detection  in  impulsive  channels  are  usually 
nonlinear,  the  analysis  of  such  problems  is  often  significantly  more  difficult  than  in  the  conventional 
Gaussian  case. 

Even  though  the  two  research  topics,  multi-user  communications  and  data  transmission  in 
impulsive  noise,  have  been  explored  thoroughly,  the  problem  of  multiple-access  communication  in 
the  presence  of  impulsive  noise  has  not  been  given  much  attention  (  see  [8-10] ).  From  the  previous 
work  on  impulsive  channels,  one  would  expect  that  in  multiple-access  communication  systems 
operating  in  impulsive  environments,  replacing  the  Gaussian  assumption  with  a  more  accurate  statis¬ 
tical  model  for  the  noise  and  then  implementing  receivers  that  take  into  account  the  characteristics  of 
that  noise  would  result  in  better  performance  than  that  achieved  by  the  usual  linear  correlator.  How¬ 
ever,  there  is  the  additional  factor  of  performance  against  the  multiple-access  interference  which  must 
be  considered  here,  and  one  might  expect  that  the  use  of  a  nonlinear  receiver  would  result  in  degrada¬ 
tion  in  this  aspect  of  performance.  Thus,  the  overall  effect  of  impulsive  and  multi-user  interferences 


on  linear  and  nonlinear  receivers  is  unclear. 

The  purpose  of  this  thesis  is  to  study  die  effect  of  the  deviations  from  the  Gaussian  channel 
noise  assumption  in  a  multi-user  environment  In  particular,  we  consider  the  adaptation  of  the  vari¬ 
ous  methods  of  performance  analysis  that  have  been  proposed  for  either  of  these  problems  to  the 
analysis  of  the  combined  problem.  In  this  study,  the  impulsive  channel  assumption  is  carried  through 
by  modeling  the  samples  of  noise  after  front-end  filtering  (  see  also  [9-10]  ).  These  samples  are 
modeled  as  independent  identically  distributed  random  variables  with  first-order  probability  distribu¬ 
tion  of  various  types  used  to  model  impulsive  noises.  The  multiple-access  capability  in  the  system  is 
achieved  by  the  direct-sequence  technique  in  which  the  spectrum  of  die  data  signal  is  spread  with  an 
assigned  "spreading  sequence"  (  see  [28-29]  and  [31-34]  ).  We  analyze  the  performance  of  the  linear 
and  nonlinear  DS/SSMA  receivers  in  the  described  environment  with  average  probability  of  bit  error 
being  the  chosen  measure  of  performance. 

Our  study  shows  that  the  linear  correlation  receiver  does  not  perform  as  well  as  the  Gaussian 
model  predicts  when  the  non-Gaussiau  noise  has  an  impulsive  nature,  even  when  the  signal-to-noise 
ratio  is  held  constant.  For  these  non-Gaussian  impulsive  channels,  hard-limiting  correlation  receivers 
are  also  examined  and  are  seen  to  perform  with  lower  error  probability  than  the  corresponding  linear 
correlation  receiver.  However,  as  one  would  expect,  the  linear  receiver  performs  better  than  the 
hard-limiting  correlation  receiver  when  the  channel  noise  is  less  structured  ( e.g.,  the  channel  noise  is 
nearly  Gaussian  ).  In  the  hard-limiting  correlation  receiver  example,  as  the  signal-to-noise  ratio 
increases,  the  hard-limiter  apparendy  is  not  as  effective  in  separating  the  two  users  as  the  linear  corre¬ 
lation  receiver  is.  However,  this  is  outweighed  by  the  improvement  against  the  impulsive  noise  as 
long  as  the  channel  noise  is  significant  Alternatively,  the  soft-limiting  correlation  receiver,  which  is 
expected  to  perform  reasonably  well  against  both  impulsive  channel  noise  and  and  multiple-access 
interference  noise,  is  studied  briefly.  Finally,  nonlinear  correlation  receivers  are  considered  in  a  gen¬ 
eral  context  and  their  asymptotic  performances  are  analyzed.  It  is  shown  that,  asymptotically  in  the 
length  of  the  spreading  sequences  and  within  some  regularity  on  the  nonlinear  element,  the  nonlinear 


4 

correlation  receiver  can  perform  with  single-user  error  probability  in  a  multi-user  environment  if  the 
spreading  sequences  are  asymptotically  ideal.  As  is  expected,  this  study  involves  evaluating 
numerous  computationally  cumbersome  expressions.  Most  of  the  effort  is  thus  concerned  with 
finding  efficient  methods  for  either  evaluating  these  expressions  exactly  or  approximating  them. 

The  organization  of  this  thesis  is  as  follows.  Chapter  2  covers  some  preliminary  results  con¬ 
cerning  the  direct-sequence  spread-spectrum  multiple-access  (  DS/SSMA  )  signaling  scheme  and  the 
model  for  the  impulsive  channel.  Chapter  3  includes  an  overview  of  some  results  on  the  performance 
of  linear  and  nonlinear  correlation  receivers  in  the  single-user  case.  Chapters  4  and  5  are  devoted  to 
the  exact  analysis  of  the  performance  of  linear  and  hard-limiting  correlation  receivers  respectively,  in 
the  impulsive  multi-user  environment  Approximations  for  the  error  probability  of  these  receivers  are 
developed  in  Chapters  6  and  7  for  systems  using  long  spreading  sequences.  To  gain  more  insight  into 
the  impulsive  multi-user  problem,  in  Chapter  8,  an  asymptotic  analysis  is  carried  out  for  both  linear 
and  nonlinear  receivers,  where  asymptotics  are  taken  as  the  spreading  sequence  lengths  increase 
without  bound.  Finally,  in  Chapter  9,  soft-limiting  correlation  receivers  as  an  alternative  for  the 


hard-limiter  are  considered. 


CHAPTER  2 


PRELIMINARIES 

In  this  chapter,  we  give  a  brief  description  of  b.  iS  VIA  signals  and  state  some  properties  of 
this  multiple-access  technique.  Then,  a  method  is  suggested  for  carrying  out  an  impulsive  channel 
assumption  in  this  environment  In  the  process,  we  introduce  some  classes  of  first-order  probability 
distribution  functions  that  are  widely  used  in  modeling  channel  noise  sources.  Among  these  non- 
Gaussian  channel  models,  the  e-mixture  of  two  Gaussian  distributions  and  the  Laplacian  noise 
models  are  the  most  interesting  in  the  context  under  study  here,  and  thus  we  will  concentrate  pri¬ 
marily  on  these  models.  Three  alternative  models,  the  Middleton  Qass  A,  the  generalized  Gaussian, 
and  the  Cauchy  noise,  are  introduced  in  Appendix  A. 

2.1.  General  DS/SSMA  System  Model 

The  model  of  the  asynchronous  binary  PSK  direct-sequence  SSMA  system  considered  here 
allows  K  users  to  share  a  channel.  The  signal  representing  the  kfll  user’s  binary  information  sequence, 
b^t),  is  a  sequence  of  unit  amplitude,  positive  and  negative,  rectangular  pulses  of  duration  T  that  can 
be  written  as 

bk(0=  £  b£>PT(t-mT),  k=l,2, ...,K,  (2.1) 


where  bm1  e  {-1,  + 1}  is  the  m*  information  symbol  of  the  k*  user  ( the  vector  of  information  sym¬ 
bols  of  the  k*  user  is  denoted  by  b(k) ) ,  and  Py(  • )  is  the  unit  rectangular  pulse  of  duration  T  defined 


(2.2) 


P-KO  = 


0  £  t  <  T, 
otherwise. 


The  data  signal  is  modulated  onto  a  phase-coded  carrier,  and  the  resulting  transmitted  signal  for  the 
k*  user  is  ( see  [26]  and  [27]) 

sk(t)  =  V2j^ak(t)bk(t)cos(0)ct  +  ek),  k=l,2,...,K,  (2.3) 

where  rck  is  the  power  utilized  by  the  k*  user,  ©c  is  the  carrier  frequency,  and  the  phase  angles 
0k,  1  5  k  £  K ,  are  not  identical  since  the  transmitters  used  in  such  systems  are  usually  not  phase  syn¬ 
chronous.  In  (2.3),  ak(  t )  is  the  code  waveform  generated  by  the  spreading  sequence  assigned  to  the 
kf*1  user.  This  spectrum-spreading  signal  can  be  written  as 


ak(t)=  £  a/k)PT.(t-jTc),  k=l,2,...,K, 
j— » 


(2.4) 


where  a^et-1,+1}  and  also  a/*^^  for  all  j  and  k  and  for  some  integer  N.  We  are  also  assuming 
that  N  is  the  least  period  of  the  sequences.  The  quantity  Tc  is  the  chip  length,  and  we  assume  that 

T=NTC  so  that  there  is  one  code  period  a^,  . ajSJ^  per  data  symbol.  For  each  signal,  sk(t), 

transmitted  by  the  k*  user  there  is  an  associated  delay  xk  for  a  given  receiver.  This  time  delay 
accounts  for  the  propagation  delay  and  the  lack  of  time  synchronism  between  transmitters.  The 
actual  received  signal  for  a  given  receiver  in  this  model  for  an  additive,  possibly  impulsive,  channel  is 
given  by 

r(t)  =  n(t)  +  £  \J2I\  bk(t-tj0)  cos  (©ct  +  4>t) , 


(2.5) 


where  n(t)  represents  the  channel  noise,  Pk  is  the  k*  signal’s  received  power,  and  where 
<J>k  =  0k  -  (OcTk ,  k  =  1 , . . . ,  K  ( see  Figure  2.1 ). 

Without  any  loss  of  generality,  the  linear  and  nonlinear  correlation  receivers  studied  here  are 
assumed  to  be  matched  to  the  first  of  the  K  signals  in  the  DS/SSMA  system;  hence,  only  relative  time 
delays  and  phase  angles  need  to  be  considered.  Therefore,  we  assume  Qj  =  0  and  Tj  =  0  in  the 
analysis  of  the  receiver  synchronized  to  the  first  user’s  signal.  (An  interesting  problem  not  con¬ 
sidered  in  this  thesis  is  the  effect  of  impulsive  noise  on  phase-locking  and  timing  acquisition  in 
DS/SSMA  systems.)  Furthermore,  there  is  no  loss  of  generality  in  assuming 
4>k  e  [  0,  2tc  )  and  xk  e  [  0,  T ),  2  £  k  £  K ,  since  we  are  concerned  only  with  time  delays  modulo  T  and 
phase  shifts  modulo  2rc. 

In  DS/SSMA  correlation  receivers  the  output  of  the  correlator  is  sampled  every  Tc  seconds. 
These  samples  are  then  passed  through  a  memoryless  nonlinearity.  In  these  nonlinear  correlation 
receivers,  the  decision  on  the  parity  of  a  data  bit  is  based  on  the  sum  of  N  samples  corresponding  to 
that  bit  taken  at  the  output  of  the  nonlinear  element  (see  Figure  2.2  ).  The  test  statistic  is  thus  written 
as 

N-l .  ...  N-l 

Yn=  Sz/)=  L  a®g(Zj),  (2.6) 

j-0  j-0 

where  g  ( • )  is  the  memoryless  nonlinearity  (i.e.,  g:  R— »R)  and  where  Zj  and  Zj(1)  denote  the  input  and 
output  of  the  nonlinearity,  respectively,  at  the  j*  sampling  instant.  The  additive  noise  in  the  channel 
is  assumed  to  have  a  symmetric  probability  distribution  and  the  average  probability  of  error  is  of 
interest;  thus,  it  is  natural  to  consider  only  nonlinearities  that  are  odd-symmetric  around  the  origin.  In 
this  general  setup  and  assuming  that  Pifb^  =  -1)  =  Pr  (b^])  =  +1)  =  1/2,  the  average  bit-error  proba¬ 
bility  using  the  statistic  of  (2.6)  can  be  be  written  as 


Pe  =  j-Pr  [YN^0|b(J1)  =  -l]  +  |pr  [YN<0|bi5)  =  +l  . 
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(2.7) 


We  will  primarily  be  interested  in  the  two  particular  cases  g(x)  =  x  ,  corresponding  to  the  linear 
correlation  receiver,  and  g{x)  =  sgn(x), 1  corresponding  to  the  hard-limiting  correlation  receiver. 
The  reasons  for  choosing  these  receivers  are  the  following.  For  AWGN  channels,  the  linear  correla¬ 
tion  receiver  provides  a  sufficient  statistic  for  the  single-user  case,  K  -  1.  Moreover,  since  linear 
correlation  receivers  are  relatively  simple  to  implement,  the  vast  majority  of  existing  direct-sequence 
spread-spectrum  systems  employ  linear  correlation  receivers  even  though  they  may  be  suboptimal  in 
a  multiple-access  environment  Here,  we  analyze  the  performance  of  this  linear  receiver  against 
impulsive  additive  channel  noise  and  multiple-access  interference  noise.  The  performance  of  the 
hard-limiting  correlation  receiver  is  also  analyzed  here  since  the  nonlinearity  is  extremely  easy  to 
implement  and  in  many  communication  problems  it  has  been  very  effective  against  impulsive  distur¬ 
bances. 

Analysis  of  the  expression  for  the  average  error  probability  given  in  (2.7)  requires  some  statisti¬ 
cal  assumptions  on  the  multiple-access  (  MA  )  interference  as  well  as  a  model  of  the  channel  noise. 
Our  statistical  assumptions  on  the  MA  interference  are  that  the  elements  b00,  Tk,  and  <j>k,  2  £  k  £  K  are 
mutually  independent  random  variables.  We  assume  that  b^0,  <  m  <  +«>,  is  a  sequence  of 

independent  data  bits  for  each  k  and  that  Pr  (b^  =  +1)  -  Pr  (b^  =  -1)  -  1/2  for  each  k  and  m.  The 
random  variables  Tk  and  <|>k,  2  £  k  £  K  are  assumed  to  be  uniformly  distributed  over  the  sets  of  their 
possible  values.  We  also  assume  that  these  variables  are  independent  of  the  channel  noise  and  of  b(1). 


■sl 


2.2.  Impulsive  Channel  Models 


Since  communication  systems  are  often  interfered  with  by  noises  other  than  the  classical  white 
Gaussian  noise,  it  is  necessary  to  consider  other  appropriate  (  and  tractable  )  noise  models.  There¬ 
fore,  in  this  section,  we  introduce  a  tractable  way  of  studying  non-Gaussian  channels  for  correlation 
receivers  such  as  that  of  Figure  2.2.  This  method  has  frequently  been  used  in  various  communication 
system  analyses. 

Our  main  assumption  concerning  the  additive,  zero  mean  channel  noise  n(t)  is  that  samples  rjj 
taken  at  the  chip  rate,  T^  after  front-end  filtering  are  independent  and  identically  distributed  (  see 
(3.1)  for  an  explicit  expression  forr|j ).  This  assumption  is  valid  when  the  noise  process  is  white  and 
Gaussian.  When  the  noise  process  is  white  but  not  Gaussian,  with  the  usual  low-pass  filtering  the 
samples  are  uncorrelated  but  not  necessarily  independent  at  the  appropriate  sampling  rate.  However, 
making  the  independence  assumption  is  still  justified  in  our  model  in  view  of  the  considerable  gain  in 
tractability  this  assumption  yields  (  see  also  [36]  ).  This  allows  us  to  study  non-Gaussian  impulsive 
noise  sources  by  modeling  the  first-order  probability  distribution  functions  of  these  independent  ran¬ 
dom  variables.  Non-Gaussian  modeling  thus  consists  of  considering  densities  that  are  either  heavier- 


tailed  or  lighter-tailed  than  the  Gaussian;  however,  our  interest  here  is  in  models  for  impulsive  non- 
Gaussian  channels,  which  correspond  to  longer-tailed  distributions. 

In  our  analysis  we  have  selected  the  first-order  distribution  of  the  random  variable  r|j  from 
classes  of  density  functions  with  applications  in  practice.  For  instance,  communication  in  the  low- 
frequency  (LF)  band  is  characterized  by  impulsive  atmospheric  noise  and  is  distinctly  non-Gaussian 
in  nature.  The  fact  that  the  receivers  discussed  here  implement  short-time  correlators  is  advantageous 
since  LF  communication  channels  have  inherently  narrow  bandwidth.  Some  of  the  models  for  noise 
sources  that  are  considered  here  have  been  used  to  describe  the  first-order  probability  distribution 
function  of  impulsive  noise  which  is  common  in  the  low-frequency  band  (  see  [8-10] ).  In  addition  to 
LF  communications,  the  usual  Gaussian  noise  assumption  is  inadequate  in  many  other 
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communication  problems.  Non-Gaussian  impulsive  models  are  necessary  for  phenomena  such  as 
atmospheric  noise  where  lightning  discharges  in  the  vicinity  of  the  receiver  cause  high  amplitude 
noise  spikes,  and  in  underwater  communication  problems  where  the  ambient  acoustical  noises  may 
include  impulses  due  to  ice  cracking  in  arctic  regions.  In  addition  to  these  natural  non-Gaussian 
noise  sources,  there  are  a  great  variety  of  man-made  impulsive  non-Gaussian  sources  such  as  automo¬ 
bile  ignition,  other  electronic  .devices,  and  heavy  electrically-powered  machinery. 

The  £-mixture  (or  £-contaminated  mixture)  is  a  commonly  used  and  highly  tractable  empirical 
model  for  impulsive  environments.  The  £-mixture  model  has  been  frequently  proposed  for  describing 
a  noise  environment  that  is  nominally  Gaussian  with  an  additive  impulsive  noise  component.  The 
model  was  first  applied  to  the  detection  of  signals  in  non-Gaussian  impulsive  noise  and  has  since 
been  used  in  various  statistical  communication  problems.  The  first-order  probability  density  function 
(  pdf )  of  this  noise  model  has  the  form 


V*)  =  (l-e)fnW  +  efi(x), 


where  £e  [0,1]  and  f„  and  ft  are  pdfs  [40].  The  nominal  density  function  fn(.)  is  usually  taken  to  be  a 
Gaussian  density  representing  background  noise.  The  impulsive  (or  contaminating)  component  of  the 
noise  is  represented  by  the  density  function  f[(.)  which  is  usually  chosen  to  be  more  heavily  tailed 
than  fn(.).  In  (2.8),  ft(.)  is  commonly  taken  to  be  Gaussian  with  large  variance,  and  this  is  the  model 
we  will  use  here.  The  ratio  of  the  variance  of  impulsive  component  to  the  variance  of  the  nominal 
one,  defined  as  y2  =  c^/Oq,  is  usually  assumed  to  be  on  the  order  of  10  and  100.  The  parameter  £  con¬ 
trols  the  contribution  of  the  impulsive  component  to  the  density  function.  In  this  thesis,  we  will  usu¬ 
ally  be  interested  in  the  effects  of  variations  in  the  shape  of  a  distribution  on  the  performance  of  vari¬ 
ous  systems  rather  than  in  the  effects  of  simple  changes  in  noise  power.  Thus,  in  this  model  we  will 
usually  vary  the  parameters  £  and  y2  with  the  total  noise  variance,  cr  =  (l-e)cn2  +  ECTf  ,  held  constant. 
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Among  other  useful  empirical  models  for  the  probability  density  function  of  r^,  the  Laplacian 


probability  density  function  is  an  example  of  a  moderately  heavy-tailed  distribution.  This  model  of 


impulsive  noise  channels  belongs  to  a  large  class  of  distributions  known  as  the  generalized  Gaussian 


class  (  see  Appendix  A  for  a  complete  description  of  the  class  ).  The  Laplacian  density  function  is 


given  by 


frn(x)_  V2o  6 


where  a2  is  the  variance. 


The  e-mixture  and  Laplacian  density  functions  introduced  above  yield  tractable  and  reasonably 


accurate  models  for  most  impulsive  noise  sources.  In  this  thesis,  we  will  examine  the  performance  of 


multi-user  communication  systems  against  both  e-mixture  and  Laplacian  noise  channels. 


SINGLE-USER  ANALYSIS  IN  IMPULSIVE  NOISE 


To  gain  insight  into  the  performance  of  multi-user  receivers  in  non-Gaussian  impulsive  noise, 
we  first  consider  communications  over  these  impulsive  noise  channels  in  the  single-user  context.  In 
particular,  we  evaluate  single-user  error  probabilities  of  linear  and  hard-limiting  correlation  receivers 
in  e-mixture  and  Laplacian  channels.  Appendix  B  contains  the  corresponding  numerical  results  for 
the  non-Gaussian  channel  examples  introduced  in  Appendix  A.  Furthermore,  in  this  chapter,  we 
introduce  our  asymptotic  analysis  by  considering  asymptotics  for  the  single-user  problem  (K=l).  In 
the  process,  we  state  some  well-known  signal-detection  results. 

Recall  from  Section  2.2  that  the  non-Gaussian  modeling  applies  to  the  samples  of  noise  after 
front-end  filtering.  These  noise  samples  can  be  expressed  as  (  see  Figure  2.2 ) 

(j+l)T. 

n'j=  J  n(t) cos(0Dct) dt  ,  j=0,l, .... N-l ,  (3.1) 

jTo 


where  the  T)'0,  ti'j,  . . . ,  T|'n_i  are  assumed  to  be  independent,  identically  distributed,  random  variables 
with  zero  means  and  variances  NqTc/4.  These  statistical  conditions  would  result  if  the  input  noise 
were  white  and  Gaussian  with  spectral  height  Nq/2  .  Keeping  these  second-order  conditions  constant 
allows  us  to  compare  the  results  obtained  here  with  those  previously  found  from  the  study  of  the  addi¬ 
tive  white  Gaussian  noise  channel. 

Within  these  assumptions,  we  can  evaluate  the  performance  in  the  single-user  case  of  linear  and 
hard-limiting  correlation  receivers  when  the  noise  samples  rj^’s  are  distributed  according  to  the 
examples  in  Section  2.2.  Note  that  the  spreading  sequence  is  irrelevant  to  the  performance  in  the 
air.gle-user  case  since  the  noise  is  symmetric  and  the  spreading  sequence  consists  only  of +  1  s  and 
- !  s.  When  K  =  1  the  average  error  probability  for  the  conventional  linear  correlation  receiver 


( see  Figure  3.1 )  can  be  computed  from  the  formula  in  [12]  as 


«• 

Pe  =  4— Jt'1  Ju-1(sin  u )  <&2(u)  du. 


(3.2) 


a  r 

where  <D2(u)  =  E{e,un°}  is  the  characteristic  function  of  the  sum  of  the  N  independent  identically 


distributed  random  variables,  and  r|0  ,r|i , . . . rij^  are  given  by  Tjj  = 


VnE^TV2 


;  j=0, 1,2,  ...N-l 


where  =  Pj  T  is  the  bit  energy  of  the  user.  Note  that  these  samples  of  noise  have  been  normalized 
N0 

to  have  variances - — ,  which  gives  the  formula  (3.2). 

2NE^> 

For  the  example  with  e-mixture  channel  noise,  the  expression  (3.2)  can  be  simplified  by  first 
noting  that  the  characteristic  function  <J>2(u)  is  given  as 


4>2(u)  =  S 
i=0 


N 


.  M  .  u2a2  P2 

(l-e)'eN-exp{ - 


(3.3) 


where  pj  is  defined  as 


-7  V'2 

N[(l-e)+er]  1 


i-KN-lJy2 


and  a  is  a  function  of  the  signal-to-noise  ratio  and  is 


defined  as  a  = 


2E^> 

No 


1/2 


Substituting  (3.3)  into  (3.2)  yields  the  error  probability  for  the  linear 


correlation  receiver  in  e-mixture  noise  as 


N 

P«  =  I 

iM) 


N 


(l-e)ieN-iQ(p,a). 


(3.4) 


Similarly,  for  Laplacian  channels,  the  error  probability  for  the  linear  correlation  receiver  is  com¬ 
puted  by  first  writing  an  expression  for  the  characteristic  function  as 


Sampler 
a  +i)Tc 


where  cr  is  the  variance  given  as - — .  The  single-user  error  probability  for  the  Laplacian  chan¬ 


nel  example  is  obtained  by  substituting  (3.5)  into  (3-2). 

To  demonstrate  the  single-user  performance  of  linear  correlation  receivers.  Figures  3.2  and  3.3 
have  been  generated  for  different  e-mixture  channels  and  the  Laplacian  channel.  The  figures  plot 
single-user  average  error  probability  versus  the  SNR  =  10  log  /  No  (  converted  to  dB  )  for  the 
linear  correlation  receiver.  In  the  e-mixture  example,  the  range  of  values  for  e  and  y2  are  those  used 
in  [40]  and  correspond  to  some  practical  examples.  Comparing  these  impulsive  non-Gaussian  chan¬ 
nel  examples  with  the  Gaussian  one,  these  curves  indicate  a  degradation  in  performance  for  both 
models  over  the  entire  range  of  interest  of  signal-to-noise  ratios,  with  fairly  large  degradation  in  some 
cases.  This  is  not  surprising  since  the  linear  correlation  receiver  is  designed  to  operate  on  the  white 
Gaussian  noise  channel.  These  figures  show  that  the  impulsive  character  of  the  channel  noise  can 
undesirably  degrade  the  performance  of  conventional  linear  receivers.  This  is  in  agreement  with  what 
was  expected,  and  what  has  been  observed  in  many  previous  studies.  One  interesting  observation 
comes  from  comparing  two  curves  in  Figure  3.3  corresponding  to  examples  with  e  -  0.01,  y2  -  100 
and  e  -  0.1,  y2  -  100.  In  these  examples  with  a  fixed  signal-to-noise  ratio,  increasing  the  amount  of 
contaminated  noise  from  e  -  0.01  to  e  -  0.1  improves  the  performance  of  the  linear  correlator.  This 
phenomenon  occurs  because  the  total  noise  variance  is  held  constant;  thus,  with  fixed  y2,  variation  in 
performance  is  not  monotonic  with  changes  in  e.  In  fact  two  channels,  one  with  e  -  0.0  and  the  other 
with  e  -  1.0,  result  in  identical  error  probabilities.  Apparently  a  breakpoint  of  the  error  probability 
versus  e  is  near  0.01. 

Typically,  optimum  and  locally  optimum  detectors  for  non-Gaussian  channels  are  obtained  by 
inserting  appropriate  nonlinearities  into  the  structure  of  the  correlator  as  illustrated  in  Figure  2.2. 


OMn 
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Figure  3.2.  Single-user  error  probability  for  the  linear  correlation  receiver  in  Gaussian,  Laplacian,  and 
e-mixture  channels,  N-31. 
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Figure  3.3.  Single-user  error  probability  for  the  linear  correlation  receiver  in  e-mixture  channels, 
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Among  the  many  possible  nonlinear  correlation  receivers,  the  hard-limiting  correlator  is  considered 
here  (  see  Figure  3.4  ).  The  hard-limiter  (sign  detector)  is  extremely  easy  to  implement  digitally  and 
introduces  virtually  no  processing  delay  since  it  only  checks  the  signs  of  the  samples.  Moreover,  in 
many  communication  problems,  hard-limiting  is  known  to  be  effective  against  impulsive  distur¬ 
bances. 


In  the  single-user  case,  it  is  straightforward  to  see  that  1/2  ( YN + N )  formed  by  the  test  statistic 
of  the  hard-limiting  receiver  is  binomially  distributed  under  either  bit  condition.  First  we  write  the 
average  bit  error  probability  for  the  hard-limiting  correlation  receiver  as 


P.-1/2  £  jpr  [yN=m|  b«J‘>=-l]+Pr  [ym«- m|  b<P=+ 1]  J  . 


(3.6) 


Since  1/2(YN+N)is  binomially  distributed,  we  have  the  following  equality 


Pr  [vN  =  m|  b^=-lU 

L  J  (N+m)/2 


p(N+m)/2  ^  j  _pj(N-m)/2  , 


(3.7) 


where  p  =  Pr[ri(1)£l]  and  where  T)(l)  is  a  typical  zero-mean  noise  sample  with  variance  -  For 

2E^n 


odd  integers  N,  substituting  (3.7)  into  (3.6),  and  noting  the  symmetry  of  the  system  the  average  bit 
error  probability  for  the  hard-limiting  correlation  receiver  in  the  single-user  case  is  written  as 


N 

Pe=  2 


N+l  I J 


N 


P'(l-P) 


(3.8) 


Figures  3.5  and  3.6  are  generated  to  show  that  hard-limiting  offers  substantial  improvement  over  the 
linear  correlation  receiver  for  the  impulsive  noise  examples  used  above  over  the  interesting  range  of 
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Figure  3.6.  Single-user  error  probability  for  the  hard-limiting  correlation 
channels,  N=31. 
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signal-to-noise  ratios  of  interest.  The  figures  plot  single-user  average  error  probability  of  the  hard- 
limiting  correlation  receiver  versus  the  SNR  for  these  noise  models.  In  Figure  3.7  the  single-user 
error  probability  of  the  linear  and  hard-limiter  correlation  receiver  is  drawn  versus  SNR  for  the  e- 
mixture  example  with  e  =  0.1  and  y2  =  100.  This  is  a  case  in  which  the  linear  correlation  receiver  suf¬ 
fered  large  degradation  in  performance  compared  to  the  Gaussian  channel  case.  The  error  probabili¬ 
ties  for  the  Gaussian  examples  are  also  included  in  the  figure  for  comparison. 

The  performance  of  linear  and  hard-limiting  correlation  receivers  in  impulsive  noise  within  the 
single-user  context  can  be  computed  via  (3.2),  (3.4),  and  (3.8).  Succeeding  chapters  of  the  thesis  are 
devoted  to  analyzing  the  performance  of  these  receivers  in  the  presence  of  both  impulsive  and  multi¬ 
user  noises. 

We  now  introduce  an  asymptotic  analysis  for  the  single-user  communication  problem  in  which 
we  investigate  the  behavior  of  direct-sequence  correlation  systems  when  the  signature-sequence 
length  used  per  data  bit  becomes  infinitely  large.  We  keep  the  bit  interval  and  the  signal  and  noise 
energies  per  bit  constant.  In  the  limit,  of  course,  these  conditions  would  usually  require  infinitely 
large  channel  bandwidth,  and  at  the  receiving  end,  samples  would  be  taken  infinitely  fast.  Also  note 
that  in  allowing,  N—H*»,  one  must  require  the  carrier  frequence  a>c  to  increase  without  bound  as  well. 
Although  these  conditions  are  not  realistic,  they  do  provide  information  about  the  limiting  behavior 
of  Pe  that  could  be  useful  in  cases  when  N  is  very  large  but  still  finite.  The  multi-user  situation  will 
be  treated  in  a  later  chapter. 

We  begin  the  analysis  with  the  usual  interpretation  of  the  binary  information  extraction  as  an 
hypothesis  testing  problem.  When  there  is  only  one  user  transmitting,  deciding  on  the  parity  of  a 
given  data  bit  is  a  choice  between  the  two  hypotheses 


S' 


Ho:vj=<l|-'V  a7T 


j  =  0,1,2 . N— 1,  (3.9) 

where  f)j  =  VN  ri'j.  Notice  that  in  (3.9),  Vj  =  VNZj  and  it  means  the  actual  received  signal  after  the 
front-end  filter  has  been  multiplied  by  VN  to  expand  data  which  otherwise  were  clustered  around  the 
origin.  This  scaling  does  not  change  the  problem  in  any  fundamental  way;  however,  the  reason  for 
this  will  become  clear  in  a  moment.  In  (3.9)  the  fy’s  that  are  distributed  according  to  examples  (2.7), 

T 

or  (2.8)  for  all  0<j£N-l  have  zero  means  and  variances  p2  =  — ,  independently  of  N.  With  this 

N-l 

view  of  the  problem,  the  test  statistic  based  on  sequence  length  N  is  YN  =  £  g(vj)  where,  for  exam- 

ie,  g(x)  equals  x  for  the  linear  correlation  receiver  and  g(x)  equals  sgn(x)  for  the  hard-limiting  corre¬ 
lation  receiver.  The  scaling  of  the  data  by  vN  allows  us  to  analyze  a  fixed  nonlinearity  g  as  N 
changes.  In  other  words,  with  T,  E^\  and  No  fixed,  the  common  variance  of  the  original  noise  sam¬ 
ples  rj'j  decreases  as  1/N .  In  order  to  consider  the  asymptotic  properties  of  a  fixed  nonlinearity,  we 
remove  this  decrease  by  scaling  rj'j  to  f)j . 

Detection  efficacy  is  a  measure  of  asymptotic  detection  performance,  which  for  the  above  prob¬ 
lem  is  defined  as  (see  [7],  and  [21-24]) 


(3.10) 


v  =  lim  i 

N-*  +  a» 


'3E[Yn  |  H,  ] 


30n 


N  Var(YN  |H0] 
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where  6n  - 


Here  E  [  YN  |  Hj  ]  is  the  expected  value  of  YN ,  given  that  the  hypothesis  H1 


is  true  and  Var[YN|Ho]  is  the  variance  of  Y^.  given  that  the  hypothesis  Hq  is  true.  One 
significance  of  detection  efficacy  is  that,  within  mild  regularities  on  the  nonlinearity  g  and  the  noise 
distribution,  it  relates  to  the  asymptotic  error  probability  via 


lim 

N— »  +  <*> 


(3.11) 


For  the  particular  nonlinearities  of  interest  here,  (3.11)  is  valid  for  all  of  the  impulsive  noise  models 


under  consideration.  The  efficacies  are  v  =  — — -  for  the  linear  correlator  and  v  =  4f^  (0)  for  the 


hard-limiting  correlation  receiver,  where  f^(.)  is  the  probability  density  function  of  the  scaled  sam¬ 
pled  noise. 

Recall  that  in  this  asymptotic  analysis  the  energy  per  data  bit  for  the  user  E^,  the  length  of  the 
data  period  T,  the  noise  power,  and  hence  the  signal-to-noise  ratio  are  kept  constant.  Table  3. 1  con¬ 
tains  the  asymptotic  average  bit-error  probabilities  for  the  linear  correlator  and  the  hard-limiting 
correlation  receiver  in  the  single-user  case  with  an  SNR  of  8  dB.  Figures  3.8  and  3.9  are  drawn  to 
exhibit  the  asymptotic  performance  of  these  receivers  over  the  interesting  range  of  signal-to-noise 
ratios  in  non-Gaussian  impulsive  noise.  These  results  indicate  considerable  improvement  in  the 
asymptotic  performance  by  using  hard-limiting  correlation  receivers  in  place  of  linear  correlators  for 
impulsive  noise  channels  for  the  entire  range  of  interest  of  signal-to-noise  ratios.  Note  that  this 
improvement  is  generally  more  pronounced  than  in  the  N  -  31  case  treated  above,  particularly  for  the 
Laplacian  case. 
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TABLE  3.1.  ASYMPTOTIC  ERROR  PROBABILITY  FOR  LINEAR  AND 
HARD-LIMITING  CORRELATION  RECEIVERS  IN  GAUSSIAN,  LAPLA- 
CIAN,  ANDe-MIXTURE  CHANNELS,  SNR=8  dB. 


DISTRIBUTIONS 

LINEAR 

HARD-LIMITER 

GAUSSIAN 

1.91  x  10-4 

2.30  x  10~3 

LAPLACIAN 

1.91  x  KT* 

2.53  x  10~7 

e-MIXTURE 

e -0.1  7* -10 

1.91  x  KT4 

1.36  x  KT* 

e  =0.01  f  =100 

1.91  x  KT4 

3.71  x  10“5 

e  =0.1  ■f  =100 

1.91  x  KT4 

8.29  x  10"' 8 
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Figure  3.9.  Asymptotic  error  probability  for  the  hard-limiting  correlation  receiver  in  e-mixture  noise 
channels. 
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CHAPTER  4 


ANALYSIS  OF  LINEAR  CORRELATION  RECEIVERS 


IN  IMPULSIVE  AND  MULTI-USER  NOISES 


In  this  chapter,  we  study  the  performance  of  linear  receivers  in  the  direct-sequence  multi-user 


environment  over  impulsive-noise  channels.  Appendix  B  contains  the  corresponding  results  for  the 


non-Gaussian  channel  examples  introduced  in  Appendix  A.  Here,  in  particular,  we  wish  to  study 


linear  digital  correlation  receivers  (see  Figure  3.1)  matched  to  the  first  of  K  users.  In  such  receivers, 


the  output  of  the  chip-matched  filter  is  sampled  every  Tc  seconds.  The  test  statistic  is  the  sum  of  N  of 


these  samples. 


N-l 

yn=  £Zj(1), 

j=0 


where  Zj(1)  is  aj(1)  times  the  sampled  output  of  the  front-end  filter.  When  the  input  signal  is  the  noisy 
multi-  user  signal  of  (2.5),  Zj(1)  is  given  by 


K  /  P  0+1)T.  /  p 

Zj0)  =  tV/D  +  2  y cos<t>k  j  a/1)  ak(t-tk)  b^t-T*)  dt  +  Y y  Tc  , 


j=0,l . N-l, 


where  double-frequency  terms  have  been  ignored,  and  is  defined  as 


=  |  n(t)  a/0  cos  coct  dt,  j=0,l,...,N-l. 
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The  noise  components  of  the  sampled  signals,  ti'^,  . . .  are  assumed  to  be  independent  and 


identically  distributed  random  variables  with  zero  means  and  variances 


N0Te 


It  is  easy  to  see  that 


Yn  given  in  (4.1)  is  not  a  sufficient  statistic  for  deciding  all  symbols  ( for  K  >  1 )  even  for  the  Gaus¬ 
sian  case.  For  this  problem,  a  more  complex  receiver  involving  sequence  detection  may  provide  a 
statistic  which  permits  a  more  nearly  optimal  decision  (see,  e.g.,  [41-42]  ).  However,  correlation 
receivers  perform  acceptably  and  are  very  simple  to  implement,  and  thus  they  have  been  used  in  most 
direct-sequence  SSMA  systems,  even  though  they  are  suboptimal  in  a  multiple-access  environment. 


To  recognize  the  contribution  of  the  non-Gaussian  noise  and  the  interference  from  other  users  in 


the  bit  error  probability,  and  since  for  the  linear  correlator  we  can  write  =  J  r(t)  ai(t)  coscoct  dt,  the 

o 


average  error  probability  of  the  linear  correlation  receiver  is  given  by  (see  [12]  for  details) 


Pe  =  ~Pr["1-T1l+Il<+1]' 


(4.4) 


A  N"1  T|'j 

where  T|,  =  £ - and  where  the  normalized  interference  term  Ij  is  given  by 


K  , -  COS<t>v 

It(b,T,<D)=£V^-^r- 

“  w  1 


b-?P  Rk.l('tk)  +  Rk.K'Cjt)  [• , 


(4.5) 


where  is  defined  as  the  ratio  of  the  powers  —  ,  for  k  -  2,  3, ....  K.  The  two  terms  Rk  j  and  Rk,] 


appearing  in  (4.5)  are  the  continuous-time  partial  cross-correlation  functions  between  the  km  and  first 
code  waveforms,  defined  by  ( see  also  [27] ) 
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I 


Rk>1(-c)  =  Jak(t^c)a1(t)dt, 
o 


(4.6) 


and 


Rk.i('C)  =  Jak(t-<)ai(t)dt, 


(4.7) 


for  0  <  T  ^  T.  For  values  of  x  in  the  range  0  £  mTc  ft$(m+  1)TC  <,  T,  the  two  cross-correlation 
functions  R^i,  and  Rk,i  can  be  written  as 


Rm(t)  =  CM(m-N)Tc  +  [Ck>1(m+l-N)-Ckii(m-N)](t-mTc) , 


(4.8) 


and 


RfciCO  =  CM(m)Tc  +  [Ck>i(m+l)-Ckii(m)](x-mTc) , 


(4.9) 


where  C^i  is  the  discrete  aperiodic  cross-correlation  function  for  the  sequences  a®  and  a(1)  defined 
by 


Ciu(m)  =  \ 


ma.(k)  ad)  0  £  m  £  N— 1 

Z*  dj  dj+m 

N"yma(k)  a(1>  ,-Nim<0 

2*  “j-m  aj 

j*° 

0  M 


(4.10) 


.  ■TO  . 


Recalling  the  statistical  assumptions  on  the  asynchronous  DS/SSMA  system,  we  note  that  I]  has 
a  symmetric  distribution,  and  since  T|,  and  Ii  are  assumed  to  be  independent,  the  average  error  proba¬ 
bility  can  be  written  as 


-re  1  Ju_1  (sin  u)  <D2  (u)  du  +  re-1  Ju-1  (sin  u)  <X>2  (u)  [1— <I»j(u)3  du, 
o  o 


(4.11) 


where  <D2  (u)  =  E  {  ehnv}  =  [  E  {  eIUT1# }  ]N  as  in  (3.2)  and  <X>i(u)  =  E{eIUl,}.  Since  Pe  is  written  con- 
veniendy  in  terms  of  characteristic  functions,  contributions  of  the  multiple-access  ( MA )  interference 
and  noise  are  distinguishable.  Note  that  the  first  two  terms  of  the  expression  for  the  multi-user  error 
probability  are  identical  to  the  right-hand  side  of  (3.2).  The  characteristic  functions  <D2  (u)  are  given 
in  the  previous  chapter  for  all  noise  models  under  consideration. 

With  regard  to  <Di  (u)  in  (4.1 1),  since  the  K-l  multiple-access  interference  terms 


b3>  Rk.KXk)  +  bj»  Ru(Tk)  >■,  2£lc£K 


(4.12) 


K 

(  i.e.,  these  form  the  total  multiple-access  interference  as  Ii  =  £Ik.i(bk ,  Tk ,  <j>t)  )  are  mutually 

te*2 

independent,  we  can  write  (  see  also  [12] ) 


4>i(u)=nj 


k=2 


2x  T 

(8rfT)-1  2  J  J  exPt  ‘“M1* .  Xk ,  4>k)  ]dxd<t> 

b-i,bo  0  0 


(4.13) 


where  £  denotes  the  sum  over  all  (bij) ,  b^* )  with  b^et+l  ,  -1}.  The  expression  in  (4.13)  can  be 

b-i,bo 


written  as 


(4.14) 


K  1 

=nir 

k-2ZXN 


Z  [f(u,m,Viu  Qk,i)  +  f(u,m,Ve][J eu)]  ; 

n*K) 


where 


J  2x  T.  r 

f(u,m,y)  =  — —  [  J  exp*  iu  [y(m+l)T  +  y(m)(Tc-^)]  Idx  d0 , 
•J**c  00  * 


(4.15) 


for  an  arbitrary  function  y.  In  (4.14)  0^  and  0^  are  the  periodic  and  odd  cross-correlation  functions 
for  the  binary  spreading  sequences,  respectively,  and  they  are  defined  as  (see  [27]  and  [33] ) 


9k,i(m)  =  Ctl(m)  +  C^ifm-N)  , 


(4.16) 


0k.i(m)  =  Ck.1(m)-Ck,1(m-N), 


(4.17) 


for  0<  ra  <  N-l  .  For  the  binary  PSK  system  considered  here,  (4.15)  reduces  to  the  following  expres¬ 
sion  which  was  obtained  in  [12]  and  [11]: 


*  2  4  y 

f(u,m,y)  =  —  Jcos[~r(y(m+l}fy(m))cos0]  sinc[— (y(m+l)-y(m))cos0  ]  d0 
tt  -  ZN  ZTtN 


(4.18) 


Typical  average  bit  error  probabilities  of  the  linear  correlation  receiver  for  the  two-user 
(K  -  2)  case  with  various  values  of  distribution  parameters  obtained  by  substituting  (4.14)  and  the 
corresponding  expression  for  <P2  in  (4.1 1)  are  contained  in  Table  4.1.  For  an  example  with  N  =  31  a 
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m 
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TABLE  4.1.  AVERAGE  ERROR  PROBABILITY  OF  LINEAR  CORRELA 
TION  RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ;  GAUSSIAN 
LAPLACIAN,  AND  e-MIXTURE  NOISE  CHANNELS,  SNR-8.0  dB. 


DISTRIBUTION 


GAUSSIAN 


LAPLACIAN 


e-MIXTURE 

0.1,10 

0.01,100 


N-31 


1.91  x  10 
4.16x10 


N=63 


3.02x10 

3.46x10 


.9 

3.06  x  1QT4 


4.02x10“* 

2.93  x  10“* 

6.76  x  10“* 

4.21  x  KT4 

1.55  x  1CT 
1.76  x  1(T 


spreading  sequence  (  also  known  as  signature  sequence  )  of  length  31  is  assigned  to  each  user 
corresponding  to  top  entries  of  Table  A.l(a)  of  [28],  These  are  auto-optimal  least-sidelobe  energy 
(AO/LSE)  phases  of  maximal-length  sequences  (m-sequence)  of  period  N.  The  results  presented  here 
are  typical  of  our  findings  for  a  variety  of  system  parameters.  For  simplicity,  die  examples  are  carried 
through  under  the  assumption  that  all  users  have  the  same  power.  All  of  these  results  are  obtained 
with  fixed  signal-to-noise  ratio,  thereby  showing  the  effects  of  the  shape  of  the  noise  density  on  per¬ 
formance.  For  the  more  interesting  example,  the  e-mixture  case,  Figure  4. 1  is  generated  to  illustrate 
the  performance  for  a  wide  range  of  values  for  SNR.  Also  shown  on  Figure  4.1  is  the  performance 
obtained  in  a  six-user  case  Gaussian  channel. 

The  principal  conclusion  here  is  that,  as  expected,  the  linear  receiver  does  not  perform  as  well 
as  the  Gaussian  model  predicts  when  the  non-Gaussian  noise  has  an  impulsive  nature  (heavy-tailed 
distribution).  Also,  the  effects  of  the  impulsive  noise  on  the  two-user  system  seems  to  be  nearly  the 
same  as  that  on  the  single-user  system.  Similar  to  the  single-user  case,  the  multi-user  performance  of 
linear  correlation  receiver  in  the  e-mixture  channel  with  e  -  0.1,  y2  -  100  is  better  than  that  in  the 
channel  with  e  =■  0.01,  y2  -  100.  Asa  final  remark  on  these  results,  in  Figure  4.1  we  see  that  for  most 
values  of  SNR  the  error  probability  for  the  linear  correlation  receiver  in  one  of  the  e-mixture  exam¬ 
ples  with  two  users  is  nearly  as  high  as  that  for  the  Gaussian  example  with  six  users.  This  indicates 
that  if  the  Gaussian  assumption  for  a  channel  is  violated  in  favor  of  this  e-mixture,  the  resulting 
degradation  in  performance  is  equivalent  to  that  caused  by  subscribing  four  additional  users  to  a 
Gaussian  channel. 

Next  we  consider  the  near-far  effects  on  the  performance  of  the  linear  DS/SSMA  correlation 
receiver  in  an  impulsive  channel.  We  assume  that  two  users  share  an  e-mixture  channel  and  that  the 

P2 

interfering  user  has  power  different  than  the  receiver  one  ( i.e.,  62,1  =  —  *  1  )•  Results  for  a  fixed 

*  1 

signal-to-noise  ratio  (  of  user  1  )  equal  to  8  dB  and  a  value  of  N=31  are  contained  in  Table  4.2  and 
Figure  4.2.  The  error  probabilities  indicate  smaller  degradation  in  performance  due  to  the  impulsive 
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Figure  4.1.  Error  probability  for  the  linear  DS/SSMA  correlation  receiver  in  e-mixture  channels, 


TABLE  4.2.  ERROR  PROBABILITIES  OF  LINEAR  CORRELATION 
RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ;  GAUSSIAN  AND 
£- MIXTURE  CHANNELS,  SNR-8.0  dB,  N-31,  K=2,  TWO  USERS  WITH 
UNEQUAL  POWERS  (  ,  =  P2IP.  SI ). 


DISTRIBUTION 


GAUSSIAN 


e-mixture 
e=0.1,  T^-IO 
e-O.OL^-lOO 
e-O.l.^-lOO 


in 


4.16  x  HT* 

7.84  x  10"4 

1.32  x  10~3 

6.76  xKT4 

1.08  x  lO-3 

1.63  x  1(T3 

3.46  xl(T3 

3.94  xl0“3 

4.49  x  lO*3 

1.31  x  ir3 

1.79  x  10”3 

2.39  x  1(T3 
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noise  component  as  the  other  user’s  power  increases  when  compared  to  the  equal  power  example. 

This  is  due  to  the  fact  that,  as  the  other  user’s  interference  becomes  more  significant,  the  channel 

noise  represents  a  smaller  fraction  of  the  total  noise.  Consequently,  the  performance  of  the  linear 

correlation  receiver  in  these  impulsive  channels  becomes  comparable  to  the  performance  of  that 

0 

receiver  in  a  Gaussian  channel  under  similar  conditions. 

Summarizing  this  chapter,  we  have  considered  the  performance  of  linear  correlation  reception 
in  a  typical  two-user  DS/SSMA  system  operating  in  impulsive  noise.  As  expected,  we  have  seen  that 
performance  in  this  situation  can  degrade  undesirably  from  that  predicted  for  the  Gaussian  noise,  par¬ 
ticularly  when  the  channel  noise  is  not  dominated  by  the  multi-user  interference. 

Aside  from  the  these  findings  it  should  be  noted  that  we  experienced  a  rapid  growth  in  the  com¬ 
putational  effort  required  to  compute  Pe  for  larger  N,  as  one  would  expect.  The  problem  of  approxi¬ 
mating  Pe  for  large  values  of  N  will  be  discussed  in  Chapter  6. 
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CHAPTER  5 


ANALYSIS  OF  HARD-LIMITING  CORRELATION  RECEIVERS 
IN  IMPULSIVE  MULTI-USER  CHANNELS 


Motivated  by  the  finding  of  the  previous  chapters,  in  this  chapter  we  study  the  performance  of 
the  hard-limiting  correlation  receiver  under  the  same  conditions  as  those  for  which  the  linear  correla¬ 
tion  receiver  was  analyzed.  We  introduce  limiting  to  diminish  the  influence  of  large  samples  (e.g., 
noise  impulses)  on  the  test  statistic.  In  Chapter  3  we  showed  that,  in  the  single-user  environment, 
hard-limiting  correlation  receivers  can  be  used  to  suppress  the  effects  of  non-Gaussian  impulsive 
noise  and  thus  to  achieve  a  lower  average  probability  of  error.  Here,  we  analyze  the  performance  of 
the  hard-limiting  correlation  receiver  in  the  presence  of  non-Gaussian  impulsive  noise  and  multi-user 
interference,  and  we  compare  these  combined  effects  with  those  on  the  linear  correlation  receiver.  As 
noted  before,  the  nonlinear  element  in  the  structure  of  the  hard-limiting  receiver  is  extremely  easy  to 
implement  digitally  and  virtually  no  processing  delay  is  introduced  since  it  only  checks  the  signs  of 
the  samples,  Zj’s  (see  Figure  3.4  ). 

The  average  bit-error  probability  for  the  hard-limiting  receiver  can  be  derived  from  the  general 
form  (2.7).  We  assume  throughout  that  the  number  of  chips  per  bit  is  odd  ( this  is  usually  the  case  as 
N  is  usually  2n-l  for  some  integer  n  )  and  so  since  sgn  (x)  e  {-1,+1},  the  test  statistic  is  an  odd 
integer,  YNe  {-N,-N+2,...,-3,-l,+l,3,...,N}. 

It  is  straightforward  to  see  that  the  test  statistic  1  /  2  (  Yn  +  N )  formed  by  the  hard-limiting 
receiver  is  multinomially  distributed  under  either  bit  condition.  From  this  it  follows  that  the  average 
Sor  probability  for  the  hard-limiting  correlation  receiver  can  be  written  as 


I  |Pr[YN  =  2m-l  |bi«=-l,x,0,bj 


+  Pr  |VN  =  -2m+l  |  b,il)  =  +l,i,<j>,b  1, 


where  x  =  fa , . . .  ,Tr  )  with  xke[0,T),  k  =  2, ....  K ,  ^  =  (<j>2 .....  <|>k  )  with  <t>ke[0,2jt ) , 
k  =  2, . . . , K ,  andb  =  [(bi?,  b^2> ) . (bif .  b^ )]  with bp)e{-l,+l }  ,k  =  2,. . .  ,K. 

We  now  consider  two  methods  for  computing  the  average  Sor  probability  (5.1)  for  the  hard- 
limiting  correlation  receiver  in  multiple-access  noise.  Results  computed  from  these  will  then  be  com¬ 
pared  with  analogous  results  for  the  linear  correlation  receiver  from  Chapter  4. 


5.1.  The  Characteristic  Function  Method 


Our  first  method  for  computation  of  Pe  for  the  hard-limiting  correlator  uses  properties  of  the 
characteristic  functions  of  discrete  random  variables.  The  characteristic  function  of  Yn  conditioned 
on  ( x ,  <t> ,  b )  and  b^1*  is  given  by 


<*>Y„  I  (U)  =  n  (eiuPjb^  +  e-iu(  1  -  p*")] 
45<r.pi  j-o 


where  Pjbi  >  =  Pr  Jzj(1)  >  0  |  b^pjj,  and  where  pi  is  used  to  denote  the  collection  of  variables 

(x,<(>,b).  Nore  that  the  probability  of  a  discrete  random  variable  taking  an  integer  value  m  is  given 
in  terms  of  its  characteristic  function  as 
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1  231 

Pr(YN  =  m)  =  — J  e_mm 0Yn(u) du; 


(5.3) 


i.e.,  it  is  the  m*  Fourier  coefficient  of  die  characteristic  function.  This  expression  can  be  simplified 
using  some  properties  of  the  characteristic  function  .  First  note  that  <DYk(u) =-<&y„(u  +  rc)  and 
Oyn(u)  =  -Oyn(ti-u)  ,  because  N  is  assumed  to  be  an  odd  integer.  Now  since  in  (5.1)  the  summation 
is  taken  over  all  the  odd  integers  m  between  1  and  N,  we  only  need  to  consider  odd  integers  m.  Equa¬ 
tion  (5.3)  can  now  be  simplified  as 


<- * 
2 


Pr(YN  =  m)  =  - 

7C 


J  e-mm  d>yN(u)  du  +  J  eiuai  dVN(u)  du 


Moreover,  for  all  odd  integers  m  it  follows  that 


(5.4) 


Pr  (yn  =  m  |  b^>=-l,pi]  =  |  jRe-e^ <X>yn(u |-1,| 


Pi) 


fdu 


Before  substituting  (5.2)  into  (5.5),  we  note  the  following 


(5.5) 


N-l 


e",um  <DYn(u)  =  n 

j-o 


r  /«  ni .  *  .  .  um  .  . 

t  cosu(l+-^-)  +  2pj  sin  u  sm(— )  ] 


-i[  sinu(l+-^-) - 2pj  sin u  cos(-^p-)  1 1, 


(5.6) 


where  Pj  =  Pjb^'  with  b<S1)  = — 1.  Substituting  (5.6)  into  (5.5)  and  after  a  few  straightforward 


*  /* ;v 


15 


intermediate  steps,  the  expression  (5.3)  can  be  written  as 


Pr  {VN  =  m  |  = — l.Pil  =  ~  j  FI  R(u,  Pj)  cos  [  £  A(u,  pj,  m)  ]  du , 

^  J  it  s_n  :_n 


o  j«0 


where 


R(u,p)  =  [1  +  4p(p-l)  sin2  u]h  , 


,,  ,  m  v  .  uni 

sinu(l  +  — -  )-2psinucos  - 

A(u,p,m)  =  tan-1 - £ - 5L 

cos  u(l  +  +  2p  sin  u  sin 

N  N 


Similarly,  we  can  write 


Pr  (VN  =  -m  |  b^^+l.pj]  =  -^-  j  ]~[  R(u,  pj) cos  [  £  A(u,^,m)]du, 

J  K  0  j-0  j-0 


(5.10] 


where  ^  =  Pr  ( Zj(1)  <  0  |  b^t}  =  +1,  pj  ]. 

At  this  point,  for  simplicity  in  the  analysis,  we  assume  that  relative  time  delays  x  among  the 
users  are  fixed.  This  analysis  can  be  applied  to  a  sufficiently  rich  set  of  delays  to  obtain  an  approxi¬ 
mation  to  the  average  bit-error  probability.  We  write  the  probabilities  pj  and  explicitly  as  functions 
x,  via 


•**  v' 1  *  v.  W’r  rv'iw 


r,ir 


***;*6*,v 


.j,*1*? 
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Pj  ( t )  =  Pr  [nj(1)  +  Ij(1)  (b,  t,  ^  1  i  x,  b,  4|] 


Vi  (x )  =  Pr  [r|/l)  + 1/‘>  (b,  x, «  <  -1  |  x,  b,  $] , 


(5.11) 


where  T|  ;(1)  =  — — — —  has  zero  mean  and  variance 

A/?-. 


'.(1) 


NoN 
2Eb(1) ' 


The  j*  sample  of  the  interference  is  written  as 


Ij(1)  (b,  x,  <j»  =  £  ““  [  B(j.ntk)  a/1)  a$*-iX'k  +  B  a/1)  a^Tc-'t'k)  ] ,  (5.12) 


fc*2 


*h“  ^  ^  ‘  ^  ■*  i-“d  -  power,  |^.  W*. *.  and 

the  functions  B(  • ,  • )  and  B  ( • ,  • )  are  defined  to  account  for  the  possible  partial  overlap  of  adjacent 
bits;  i.e.. 


BO.mt)  = 


bi^  0£j:Smk 
btf*  mk+l  £  j  £  N-l 


and 


B  (j.nik)  =  < 


blf*  0£jSmk-l 
b^  mk£jSN-l 


(5.13) 
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Substituting  (5.7)  and  (5.10)  into  (5.1),  the  average  bit-error  probability  conditioned  on  the 
delays  is  written  as 


i  (N+l)/2* 


Recall  that  the  tj/^’s  arc  i.i.d.  with  a  symmetric  density  function.  Thus 


N-l 


N-l 


n  R(u,pj(D)cos[  2  A(u,Pj(r),2m-l)] 


o  U-o 


j-0 


du 


d<j> . 


(5.14) 


Pr[T|/1)  S  1  -  Ij(1)]  =  Prfo/”  <  -1  + 1/1)],  and  T,  b)]  =  0. 


Now  since  (5.7)  and  (5.10)  are  functions  of  pj(t)  and  ^(x),  respectively,  and  with  the  above  properties 
we  can  express  P„(x)  as  in  (5.14). 

In  the  case  that  only  two  users  share  the  channel  (K=2),  (5.14)  can  be  simplified  significantly  by 
noting  that  the  absolute  value  of  Ij(1)  takes  only  two  forms: 


l/1>(b,T,<t>)  = 


f  ±  cos  <j>2  (l-2d2) 

±  Videos  <(>2, 


(5.15) 


where  d2  =  t'j/T,..  Now  defining  jj,  i=l,2,3,4  as  the  number  of  times  in  N  tries  that  Ij(1)  takes  one  of 
the  four  forms  (e.g.,  Ij(l)  =  cos  4>2( 1— ^d2))  enables  us  to  simplify  (5.14)  as 


Pe(I) 


4  .  (N+iy2  4 

IIRCu.pWy1  Z  cos  [  Z  JiA(u,  p(i),2  m  - 1 )  ] 

i«l  real  i»l 


du 


d<(> , 


(5.16) 


where  p(i)  =  Pr[rjj<1)  >  1  -  when  Ij(1)  takes  the  i*  form  i=l, 2,3,4.  Note  that  j; ,  which  is  a  function 


of  the  data  bits  of  the  other  user,  also  depends  on  the  signature  sequences  of  the  two  users.  An 
expression  for  Pe(x)  similar  to  (5.16)  can  be  obtained  when  more  than  two  users  share  the  channel. 
However,  illustration  of  the  final  expression  would  be  somewhat  difficult 

5.2.  The  Combinatorial  Method 

We  now  derive  an  alternative  expression  for  die  average  bit-error  probability  from  the  original 
definition  (5.1)  when  only  two  users  share  the  channel.  Recall  that  the  objective  is  to  evaluate 
Pt[Yn  =  m]  for  some  odd  m,  where  YN  is  the  sum  of  N  random  variables  Zj(1)  e  {-1,+1}.  First  we 
rewrite  the  four  probabilities 

P(  1)  =  1  -  V^cosfc] , 

P(2)  =  PH<)*1  +^i‘co*d. 

P(3)  =  Prft/”  >  1  -  cos4>2(l-2d2)j , 

and  P(4)  =  Pr[r)j(1)  S  1  +  -\/e2.i  cos<{>2(  l-2d2)J 

which  are  conditioned  on  x  and  The  test  statistic  YN  can  take  on  any  odd  integral  value  between 
-N  and  +N.  YN  equals  an  odd  integer  m  only  if  the  number  of  times  in  N  tries  comes  out  non¬ 
negative  is  exactly  (N+m)/2.  The  probability  of  each  Zj(1)  being  non-negative  is  equal  to  P(i)  for 
some  i,  i=l,2,3,4,  depending  on  the  bits  of  the  other  user  and  the  signature  sequences.  Recalling  the 
definition  forjj,  i=l,2,3,4  we  can  write 
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Pr[YN  =  m|b(J1>  =  -l,pi]  =  Z;ZE 

*■  J  it.  n«  n. 


04  0)  0] 


N+ra 


*  ~  Qa  •  n  j  —  04 


Ji  jj  j« 

0j  0]  Ot 


Mho  N+m 

— - Oa  —  lb  —  04  ji - - —  +  m  +  0)  +  n« 

PO)  2  [1-P(1)1  2  p(2)^[l-p(2)]Jj-nj 


po^ti-po)}**-”*  p(4r  [i — p(4)j**— 114 , 


where 


=  In  (5.14)  the  summations  are  in  the  ranges 


max 


N+m 


-j3“j2-ji.O 


<;  n4  S  min 


N+m 


>J4 


max 


N+m 


-n4-j2-ji.O 


5  n3  5  min 


N+m 


-  »Ud3 


and 


max 


N+m 


-114-03-^,0 


<n2  S  min 


N+m 


~  °4  ~  °3  J2 


Now  substituting  (5. 17)  into  (5.1),  we  can  write 


1  (N+iy2  2  r  n 

pe(D  =  T~I  s  JPr  YN  =  2m-l  |  b^=-l,T,b,  <j> 

b  m-1  0  L  “j 


(5.17) 


(5.18) 


As  with  the  final  simplification  (5.16)  in  Section  5.1,  the  analytical  result  obtained  here  can  be  gen¬ 
eralized  to  the  case  when  more  than  two  users  share  the  channel,  albeit  it  is  somewhat  messy  to  do  so. 
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5J.  Numerical  Results 


The  numerical  examples  presented  here  are  aimed  at  showing  that,  for  fixed  SNR,  improvement 
in  DS/SSMA  performance  can  be  obtained  by  using  hard-limiting  correlation  receivers  in  place  of 
linear  correlators  when  the  tails  of  the  noise  distribution  are  heavy.  Thus,  we  compute  average  bit¬ 
error  probabilities  given  a  relative  time  delay  between  users  for  impulsive  noise  sources  interfering 
with  the  two-user  binary  PSK  direct-sequence  SSMA  system  considered  in  the  previous  chapter. 

This  analysis  is  first  carried  out  with  fixed  signal-to-noise  ratio  and  equal  signal  power  assump¬ 
tion  for  all  the  users.  With  the  variance  of  the  random  variable  held  constant,  changes  in  error 
probability  are  analyzed  for  different  noise  distributions.  Tables  5.1  and  5.2  contain  the  average  bit¬ 
error  probabilities  for  Gaussian,  Laplacian,  and  E-mixture  examples  with  a  typical  delay.  Results 
indicate  significant  improvement  in  performance  by  using  hard-limiting  correlation  receivers  in  place 
of  linear  correlators  for  more  impulsive  noise  channels.  Moreover,  this  improvement  becomes  more 
visible  as  the  length  of  the  signature  sequences  used  by  the  channel  subscribers  increases.  Figures  5.1 
and  5.2  support  this  conclusion  over  a  range  of  SNR  values  for  the  e-mixture  example.  An  interesting 
conclusion  here  is  that,  unlike  the  linear  correlation  receiver,  degradation  due  to  the  other  user  is  no 
longer  uniform  in  the  range  of  signal-to-noise  values.  Thus,  as  the  SNR  increases,  the  hard-limiter 
apparently  is  not  as  effective  in  separating  the  two  users  as  the  linear  correlation  receiver  is.  How¬ 
ever,  this  deficiency  is  outweighed  by  the  improvement  against  the  impulsive  noise  as  long  as  the 
channel  noise  is  significant  ( i.e.,  the  SNR  is  moderate  ). 

Finally,  we  consider  the  near-far  effects  in  the  performance  of  hard-limiting  DS/SSMA  correla¬ 
tion  receivers.  In  this  context,  we  vary  the  signal  power  of  the  second  user  with  respect  to  the  first 
user’s  power  level  and  then  compute  the  error  probability.  For  Gaussian,  Laplacian,  and  e-mixture 
channels.  Tables  5.3  and  5.4  are  generated  for  values  of  €2,i  less  than  unity  and  greater  than  unity, 
respectively.  Additionally,  Figure  5.3  is  drawn  to  illustrate  the  performance  variations  over  a  range 
of  values  for  e2.i-  Comparing  these  results  with  the  analogous  ones  for  the  linear  correlator  of 
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Figure  5.1.  Error  probability  of  the  hard-limiting  correlation  receiver  in  e-mixture  channels,  N=31, 
two  users  with  equal  power  and  typical  time  delays. 
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Figure  5.2.  Error  probability  of  linear  and  hard-limiting  correlation  receives  in  e-mixture  channels, 
N-31,  two  users  with  equal  power  and  typical  time  delays. 
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TABLE  5.4.  ERROR  PROBABILITY  OF  HARD-LIMITING  CORRELATION 
RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ;  GAUSSIAN  AND 
IMPULSIVE  CHANNELS,  SNR=8.0  dB,  N-31,  K=2,  TWO  USERS  WITH 
UNEQUAL  POWERS  AND  TYPICAL  TIME  DELAYS  ( e2 1  =  P2/Pi  >1 ). 


DISTRIBUTION 

€■>  i  =  3.01 

WBfBSSHUM 

GAUSSIAN 

4.63  x  lO-3 

8.30  x  10-3 

1.32  xlO"2 

LAPLACIAN 

8.67  x  KT4 

4.05  x  ir3 

9.58  x  ir3 

e-MIXTURE 

e=0.01,f=100 

3.90  x  1(T* 

1.74  x  10r3 

4.93  x  ir3 

e-0.1,^-100 

2.59  x  KT4 

4.21  x  10"4 

3.85  x  10-3 

Chapter  4  we  see  that,  as  the  other  user’s  relative  power  increases,  there  is  a  point  beyond  which  the 
linear  correlator  is  superior  to  the  hard-limiting  correlator.  This  again  is  due  to  the  fact  that  the  linear 
correlator  appears  to  be  more  effective  than  the  hard-limiter  is  in  dealing  with  the  MA  interference. 
So  when  the  MA  interference  dominates  the  impulsive  interference,  the  linear  correlator  is  better. 
However,  this  crossover  point  is  higher  for  more  impulsive  models  and  in  most  cases  occurs  beyond 
the  desirable  range  of  error  probabilities.  Thus,  we  can  conclude  that  the  hard-limiter  correlator  is 
still  desirable  in  a  two-user  impulsive  channel  provided  the  second  user  is  not  dominant 

Summarizing  this  chapter,  we  have  considered  die  performance  of  hard-limiting  correlation 
receiver  in  a  typical  two-user  DS/SSMA  system  operating  in  impulsive  noise.  We  have  seen  that  in 
most  cases,  the  hard-limiter  outperforms  the  linear  correlator,  particularly  when  the  MA  interference 
is  dominated  by  the  impulsive  channel  noise.  Similar  to  the  analysis  for  the  linear  receiver,  we 
experienced  a  rapid  growth  in  the  computational  effort  required  to  compute  Pe  for  larger  N.  The 
problem  of  approximating  Pe  for  large  values  of  N  will  be  discussed  in  Chapter  7. 


CHAPTER  6 


LINEAR  CORRELATION  RECEIVERS  IN  IMPULSIVE  CHANNELS: 

APPROXIMATIONS 

In  earlier  chapters,  expressions  for  average  bit-error  probability  were  obtained  for  both  linear 
and  hard-limiting  correlation  receivers  operating  in  the  presence  of  multi-user  noise  and  non- 
Gaussian  impulsive  noise.  The  computational  effort  in  finding  P0  grows  rapidly  as  N,  the  number  of 
chips  per  data  bit,  and  K,  the  number  of  users  sharing  the  channel,  increase.  In  fact,  in  our  previous 
analysis  we  have  been  limited  by  computational  considerations  to  considering  exclusively  the  case 
K-2.  In  this  chapter,  we  consider  further  the  performance  of  linear  correlation  receivers  in  the 
prescribed  environment.  Maintaining  the  same  system  model,  we  introduce  an  approximation,  an 
upper  bound  and  a  lower  bound  for  the  average  probability  of  error  that  allow  for  easier  computation 
with  large  N.  The  approximation  to  the  performance  of  the  linear  correlation  receiver  is  based  on  a 
Taylor  series  expansion  of  the  average  error  probability.  Similar  ideas  have  been  used  in  [1 1,12]  and 
[43].  We  will  also  show  that  this  approximation  is  exact  asymptotically  in  N.  Furthermore,  we  con¬ 
sider  the  moment-space  bounding  technique  proposed  in  [44-45],  and  [17]  to  provide  upper  and  lower 
bounds  on  the  average  probability  of  error  in  this  case. 

6.1.  Taylor  Series  Approximation 

In  this  section,  we  obtain  an  approximation  for  the  error  probability  of  the  linear  correlation 
receiver  by  truncating  a  Taylor  series  expansion.  It  is  easy  to  show  that  with  the  test  statistic  for  the 
linear  correlator,  the  average  bit-error  probability  for  the  receiver  can  be  written  as  (  see  also  [12] ) 


(6.1) 


Pe  =  Et,^,b‘j 


where  TN<a)  =  Pr  [  £n"j  >a  ]  with  zero  mean  samples  of  noise  T|"j  normalized  to  have  variance  — . 

j-o  N 

The  multiple-access  interference  in  (6.1)  is  defined  as 


K  _  COS0V 

Ii(  b ,  x .  <D)  =  £  -=-=■  h#  R^iftk)  +  btf»  Rk.Ktk)  • , 

k-2  1 


where  ej^j  the  ratio  of  the  powers  — ,  k  =  2,3,.. .,K,  and  the  continuous-time  partial  cross-correlation 

pl 

functions  R^  and  R^j  are  given  in  (4.8)  and  (4.9),  respectively.  Following  [11,12]  and  [43]  we 
expand  the  function  'Fj*  in  (6.1)  in  Taylor  series  about  a  function  of  the  single-user  SNR, 

A 

a~y  ~n —  '  ^en  “  we  ****  expectation  of  each  term  ( assuming  this  interchange  is  permissi¬ 
ble  )  we  will  have  the  following: 


Pe  =  ‘FN(a)+  £  lM4^m>(a)  , 


where  vF1sf’)  is  the  mlh  derivative  of  4^  .  Note  that,  if  'Fn  is  analytic  the  series  in  (6.3)  converges 
uniformly  in  a  because  the  interference  is  bounded.  In  (6.3)  the  first  term  of  the  series  is  the  error 
probability  in  the  absence  of  multiple-access  noise.  The  rest  of  the  series  is  a  function  of  signal-to- 
noise  ratio,  the  shape  of  the  noise  distribution,  and  the  even  moments  of  the  multiple-access  noise. 


The  moments  E  [  if®  ]  can  be  evaluated  by  means  of  a  recursion  (see  [11,12]  and  [43] ).  We 
first  define  the  random  variables 


,.1  =  — bi?>  Rq.tC'tq)  +  bfft  Rq.K'tq)  • ,  q  =  2, 3 . K. 

i  « 


Let  Iq,i  be  the  partial  sum  of  the  interference  of  users  2,  3, ... ,  and  q.  This  quantity  can  be  written 


Iq.i=  X  Wk,l  cos  <$>k  =  Iq — l  +  Wql  cos  <t>q  ,  q  =  3,4 . K, 

k-2 


with  hi  -  W^i  cos  <t>2-  Note  that  Ij  =  IK>i,  and  since  <|>q  and  Wq  i  are  independent  we  can  write 


E[  I2?  ]  =  £  [2(m-i)l  2~2^  E[  I2*,.,  ]  E[  ],q  =  3 . K. 

v  21  J  v.  m-i  J 


(Notice  that  due  to  symmetry  of  the  multiple-access  interference,  the  odd  moments  of  Iq,i  are  zero.) 
For  computing  the  moments  of  the  multiple-access  interference  we  must  first  evaluate  E  [  W^1  ] .  To 
compute  these  moments,  we  use  a  modified  version  of  the  characteristic  function  of  Ij(  T ,  <j> ,  b ) 
given  in  (4.14).  This  yields 


i  N-l  _  _ 

E  [  Wqj ®  ]  =  — — —  £  [  fm(n>  Veq,l  ®q,l)  ^m(n*  q.l  ®q,l)  1  » 


where 
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y(n-H)  *»'  y(n)  ,2nHl 
f  /  ^  N  J  1  N  J 

m(D,y)  _ _ [  y(n-H)  -  y(n)  1 


i 

i 

i 


(6.8) 


fory(n)  *y(n+l)  and  fm(n,  y)  =  [-^]  if  y(n)  =  y(n+l) . 

Referring  to  expression  (6.3),  our  proposition  is  that  if  the  signature  sequences  are  chosen 
correctly,  then  when  N  is  large  the  effect  of  higher  degree  moments  of  Ij(b ,  x ,  <j>)  in  (6.3)  becomes 
negligible  (note  that  in  most  cases  we  can  safely  assume  that  I]  is  bounded  from  above  by  1).  There¬ 
fore,  only  a  few  terms  of  the  infinite  series  should  be  needed  to  obtain  a  small  truncation  error  for  Pe, 
and  thus  in  order  to  evaluate  the  estimate  for  Pj,  we  primarily  need  to  study  the  distribution  function 
'Fn  and  its  lower-order  derivatives. 

Recall  that  our  first-order  non-Gaussian  channel  model  describes  the  noise  samples  q".  Let  us 
start  with  the  Gaussian  distribution;  of  course,  the  sum  of  N  independent  zero  mean  Gaussian  random 

variables  each  with  variance  tt  is  a  Gaussian  random  variable  with  zero  mean  and  unit  variance. 

N 

That  is, 

eo 

'PN(a)  =  Q(a)  =  J  -  dx  .  (6.9) 

We  can  write  the  nth  derivative  of  Q  in  terms  of  the  Hermite  polynomials  Hm(x)  which  are  defined 
by  ( see  [25] )  Ho(x)  =  l,Hi(x)  =  2x,  and 

HnHl(x)  =  2xHtn(x)-2mHm_,(x) ;  m^l.  (6.10) 


The  resulting  expression  for  Pe  is 


H2l*-i(:gOE[I,2m]  . 


(6.11) 


Pe  =  Q(a)  +  ^-I 

yK  m-l 


a 

.  V2 


The  convergence  of  this  infinite  series  and  an  upper  bound  on  the  truncation  error  that  results  when 
the  series  is  approximated  by  the  sum  of  a  finite  number  of  terms  are  considered  in  [12],  [43],  and 
also  [11]. 

For  the  e-mixture  example,  recall  from  Chapter  3  that  the  probability  distribution  function  of 
the  sum  of  N  independent  random  variables,  each  having  a  density  given  in  (2.7),  can  be  written  as 

vFN(a)  =  £  [n]  (l-e)'eN~i  Q(fta)  ,  (6.12) 

i-o  u  J 


where  Pi  =  ^  ’  The  resul6ng  expression  for  Fe  for  the  e-mixture  case  is  thus 


.iM.  2m 


N  n  e  2  - 

P.  =  S  [s]  (l-ev  e-  Q(fta)  + V  2  -JST 


„  2  -  tHi  «-/  rt 


(6-13) 


and  the  truncation  properties  will  be  similar  to  those  for  the  Gaussian  case  discussed  in  [1 1]  and  [12]. 

For  this  case,  by  using  only  two  terms  of  the  infinite  series  in  (6.11)  and  (6.13)  we  obtain  an 
estimate  of  Pe  for  the  linear  correlator.  Typical  results  showing  the  accuracy  of  this  approximation 
are  depicted  in  Figure  6.1.  These  approximations  are  computed  for  Gaussian  and  e-mixture  examples 
with  fixed  signal-to-noise  ratio  (  SNR  -  8  dB  )  for  some  values  of  spreading  sequence  lengths  N.  As 
we  anticipated,  the  estimate  results  in  values  for  the  error  probability  closer  to  the  exact  ones  as  N 
grows.  For  the  N-31,  K«6  Gaussian  case,  this  method  is  less  accurate  than  for  K-2.  However,  for 


highly  impulsive  examples,  the  approximation  is  very  accurate  even  for  small  N.  This  is  most  likely 
due  to  the  small  values  of  the  higher-order  derivatives  of  heavily  tailed  distributions  at  a  Consider¬ 
ing  the  above  finding,  the  Taylor  series  approximation  appears  to  be  very  appropriate  for  estimating 
the  error  probability  for  the  impulsive  non-Gaussian  channels.  Using  this  technique,  we  approximate 
the  error  probability  of  the  linear  correlation  receiver  in  e-mixture  channels  for  three-and  six-user  sys¬ 
tems.  The  results  are  contained  in  Table  6.1  with  signal-to-noise  ratio  value  of  8  dB  and  signature- 
sequence  lengths  N-3 1  and  N-63.  As  we  learned  from  the  above  discussion,  the  estimates  in  Table 
6. 1  must  be  more  accurate  for  the  cases  with  larger  values  of  N  and  more  impulsive  noise  channels. 
We  also  see  that  when  the  multi-user  interference  becomes  more  significant  (  e.g.,  K  -  6  ),  the  chan¬ 
nel  noise  represents  a  smaller  fraction  of  the  total  noise.  Consequently,  the  performance  of  the  linear 
correlation  receiver  in  these  impulsive  channels  becomes  comparable  to  the  performance  of  that 
receiver  in  a  Gaussian  channel  under  similar  conditions. 

6.2.  Asymptotic  Analysis 

Equation  (6.3)  yields  a  useful  approximation  to  the  average  bit-error  probability  for  linear  corre¬ 
lation  receivers  in  impulsive  noise  for  large  N.  Now  we  consider  the  limitation  of  the  system  as  N 
grows  without  bound.  Recall  from  the  single-user  analysis  of  Chapter  3  that  in  considering  asymptot¬ 
ics  we  wish  to  keep  the  energy  per  data  bit,  the  length  of  the  data  period  T,  the  noise  power,  and  con¬ 
sequently,  the  signal-to-noise  ratio,  constant  as  N  changes.  The  statistic  for  the  binary  decision  for 
the  linear  correlator  is 


where  r|'j(1)’s  ,j  =  0, 1, . . .  ,N-1  are  the  noise  components  of  the  sampled  signal  with  zero  means  and 
N0TC 

variances  — — — .  The  I'j(  )’s  are  samples  of  the  multiple-access  interference  and  are  defined  as 
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TABLE  6.1.  TAYLOR  SERIES  APPROXIMATION  FOR  ERROR  PROBA¬ 
BILITIES  OF  LINEAR  CORRELATION  RECEIVERS  IN  THE  BINARY  PSK 
DS/SSMA  SYSTEM  ;  GAUSSIAN  AND  e-MIXTURE  CHANNEL,  SNR-8.0 
dB. 
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where  Ij(1)  is  given  in  (5.12).  Referring  to  (6.14),  conditioned  on  the  multiple-access  parameters 
Pi  =  (t,<t>,b) ,  the  test  statistic  Yn  is  seen  to  be  the  sum  of  N  independent  random  variables.  Thus, 
the  asymptotic  error  probability  of  the  linear  correlation  receiver  can  be  obtained  by  applying  the  cen¬ 
tral  limit  theorem  to  the  sum,  as  follows: 

With  the  test  statistic  as  given  in  (6.14),  the  asymptotic  error  probability  is  given  as 


lim 

N  — »+*• 


i  'APr(YN>0  |b^1)=-l,pI)  +  1/iPr(YN<0  |  b^+l  ,Pl)  k 


(6.15) 


This  expression  can  be  simplified  using  the  assumed  symmetry  as 


Pe=  lim  E^b- 

N— »+«• 


Pr(  Yn  <  0  |  b^l)=+l  ,pi) 


(6.16) 


Now  the  term  inside  the  expectation  on  the  right-hand  side  of  (6. 16)  can  be  written  as 
Pr(YN<0  |b(J,>=+l,p1)  = 


Pr 


YN-E{YNlb^=+l,Pl>  ^  ~ E { Yn  |  bj^=+l , pi }  ^ 

VVar(YN  |b^=+l,Pl)  <  VVar(YN|b^  =+l,Pl)  +  ’Pl 


(6.17) 


Note  that  the  conditional  mean  of  the  test  statistic  YN  is  given  by 


E  | YN  |  b^Wfl , p,  »=  a/y"tc  f  g(  Ija)  +D3, 


(6.18a) 


N-l 


where  Ij(1)  is  defied  in  (5.12)  and  relates  to  the  total  multiple-access  interference  as  Ij  =  1/N  £ 

i-o 

Similarly,  the  conditional  variance  is  given  by 


f  ,  )  NNq Tc 

Var  |yn  |  btfWl  ,Pl  j  = - j-A  (6.18b) 

Thus,  assuming  that  lim  Ii  ( x ,  4> ,  b )  exists  almost  surely,  the  central  limit  theorem  implies  that 
N-*+«»  —  “ 

N1inJ«Pr  [Yn<0  I  b^1)=+1*Pi]  =Q  (»•*•).  (6.19) 


Therefore,  since  Q  £  1 ,  the  bounded  convergence  theorem  implies  that  the  asymptotic  multi-user 
error  probability  of  the  linear  correlation  receiver  can  be  written  as 

Pe  =  ET>(il>b  Jq  fa[l+  lim  Ii(b  ,  t  ,  d> )  ] 

-  -  -  I  1  N-*+—  — 


where  Ii(b  ,  x ,  <j>)  and  a  are  as  defined  in  (6.2)  and  (6.3). 

From  (6.20)  we  see  that  single-user  performance  in  a  multi-user  environment  is  achieved 
asymptotically  by  the  linear  correlation  receiver  if  the  following  asymptotic  condition  is  satisfied: 

lim  It(b  ,  x  ,  <S>)  =  0  almost  surely. 

N-»+«»  —  — 


(6.21) 


Furthermore,  in  Chapter  8,  we  will  show  that  (6.21)  is  satisfied  if 

lim  max-^-|0ki(m)|  =0and  lim  max -77 |0i- i(m)|  =0  ,  k  =  2,3,...,K.  (6.22) 

N  -♦+••  m  N  ’  N  — »+«■  m  N 

The  existence  of  binary  sequences  satisfying  these  conditions  follows  by  applying  (6.22)  to  infinite 
sequences  proposed  by  Schneider  and  Orr  [34].  Schneider  and  Orr  show  the  existence  of  a  sequence 
set  whose  autocorrelation  and  cross-correlation  functions  approach  ideal  behavior  (6.21)  with  increas¬ 
ing  N.  Moreover,  they  prove  that  the  cardinality  of  such  sets  of  infinitely  long  sequences  grows  with 
N.  Thus,  at  least  one  large  set  of  binary  sequences  having  asymptotically  ideal  correlation  properties 
exists.  However,  by  utilizing  the  sequences  in  [34],  the  convergence  of  (6.21)  to  zero  is  very  slow. 
This  induces  a  limit  on  the  number  of  subscribers  to  the  channel  with  single-user  performance.  From 
the  discussion  in  [34],  it  follows  that  one  can  find  at  least  N1/4  sequences  such  that  (6.21)  is  satisfied. 

Figures  3.8  and  3.9  in  Chapter  3  contain  the  asymptotic  bit-error  probability  for  the  linear  corre¬ 
lation  receiver  under  the  assumption  that  (6.21)  is  satisfied,  since  the  multi-user  impulsive  noise  per¬ 
formance  in  this  case  is  the  same  as  in  the  Gaussian  single-user  case.  Note  that  the  above  asymptotic 
results  also  indicate  that  the  approximate  error  probability  obtained  from  truncating  the  infinite  sum 
in  (6.3)  is  asymptotically  exact  for  appropriately  chosen  sequence  sets. 

6.3.  Moment-Space  Bounds 

In  this  section,  we  apply  the  moment-space  bounding  technique  to  bound  the  average  bit-error 
probability  of  the  linear  correlation  receiver  in  non-Gaussian  impulsive  noise.  This  method,  which 
was  developed  by  Yao  in  [44]  and  [45],  is  based  on  a  result  from  the  theory  of  games.  To  present  the 
application  of  the  bound  to  our  problem,  we  first  state  a  simple  version  of  an  isomorphism  theorem  in 
moment  spaces  which  provides  relations  among  arbitrary  moments  of  a  random  variable. 


Note  from  (6.1)  that  our  main  objective  here  is  to  evaluate  the  mean  of  a  function  of  the  random 
variable  Ii(  t ,  <|> ,  b  ) .  Although  a  closed-form  expression  for  the  probability  distribution  function  Fj 
of  this  random  variable  is  not  available,  we  know  that  its  support  is  confined  to  a  finite  closed  interval 
E  =  [-D,D].  Referring  to  (6.2)  and  from  [28]  we  have  the  result  that  for  rectangular  chip 
waveforms 


f 

1  K 

D  =  —  Yj  max-  IC^m)!  +  |Ci^1(m  —  N) |  :  0<m<N  >■ . 

^  k=2  m  1 


(6.23) 


Now  let  Xi  and  Xz  be  a  pair  of  continuous  functions  defined  on  the  interval  E.  The  generalized 
moments  of  the  random  variable  Ii  induced  by  the  functions  Xi(Ii)  are 


mi  -  J  Xi(Ii)  dFj  —  E* ,  0 .  b 


i  =  1,  2  . 


(6.24) 


For  a  given  pair  of  continuous  functions  (  Xi  >  X2 )  >  we  denote  the  moment  space  M  as 


Af  = 


|m  =  (m1,m2)sR2  |  mj  =  J Xi( I )  dFi , i=l .2  ,F1eF(H)|, 


(6.25) 


where  P  (  E  )  is  the  set  of  probability  distribution  functions  defined  on  E  =  [  -D  ,  D  ].  It  is  easy  to  see 
that  M  is  a  closed,  bounded,  and  convex  set  in  R2.  The  isomorphism  theorem  states  that  if  we  plot 
X:(y)  versus  X:(y)  >n  R2  for  y  ranging  over  E  ,  then  the  convex  hull  H  of  the  resulting  plot  is  the 
moment-space  XI  ( i.e.  H  -  M ). 


The  application  of  this  result  to  bounding  the  error  probability  of  the  linear  correlation  receiver 
;s  given  as  follows.  Let  the  function  X;  be  defined  by  the  expression  inside  the  curly  brackets  in 
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X2(y)  =  ^N[ct(l+y)]  .  (6.26) 

I 

Now  let  Xi(y)  be  some  continuous  function  of  y  whose  generalized  moment  mj  may  be  readily 
evaluated.  If  this  function  is  carefully  chosen  so  that  the  resulting  curve  ( Xi(y) ,  foCy) )  has  a  thin 
convex  hull  over  the  range  S  ,  then  we  can  use  our  knowledge  of  mx  to  obtain  tight  upper  and  lower 
bounds  on  m2  ,  since  H  •M  and,  for  this  %2> 


m2  =  /X2(Ii)  dFj,. 


(6.27) 


We  would  like  to  choose  Xi  so  that  the  convex  hull  of  the  curve  of  fo  versus  Xi  is  thin.  How¬ 
ever,  in  choosing  Xi ,  we  must  trade  off  the  thickness  of  H  with  the  evaluability  of  mj.  A  useful  fam¬ 
ily  of  candidate  functions  for  Xi  consists  of  exponentials  of  the  form 

Xi(y)  =  exp{h[l+y]}  ,  (6.28) 


where  h  is  a  real  parameter  which  can  be  chosen  to  tighten  the  bounds.  For  each  choice  of  the  kernel 
we  need  to  evaluate  the  corresponding  moment.  Referring  to  (6.2),  we  note  that 


mi=ehn 

k-2 


exP  [  h  Ijc,i(  bk ,  Tk ,  <pjc )  ] 


(6.29) 


where  is  the  multiple-access  interference  corresponding  to  the  k*  user.  The  moment-generating 
function  of  the  random  variable  I|t,i  appearing  in  (6.29)  can  oe  written  as 


E  jexp  [  h  Ik,^  t>k , ■ xk ,  <|>k )  ] ||= =  J  I0(  hy[e^J [  -^■0ic.i(m)  +  -^j-e^m+l)  ] ) 

+  Io(hVeu['^®u(ra)  +  ^$k,i(m+l)])  dt  ,  (6.30) 

where  Iq  is  the  modified  Bessel  function  of  order  zero.  The  moment  of  Xi  (I|)  is  obtained  by  substi¬ 
tuting  (6.30)  into  (6.28). 

Having  chosen  Xi  to  be  of  the  form  (6.28),  we  use  a  graphical  technique  to  evaluate  the  upper 
and  lower  bounds  on  the  error  probability,  in  which  we  plot  Xi(y)  versus  X2(y)  over  die  range  S  and 
complete  the  convex  hull  H  graphically.  Since  we  can  compute  mi  via  (6.2S)  the  upper  and  lower 
bounds  on  the  error  probability  can  easily  be  read  from  the  graph  by  drawing  a  line  perpendicular  to 
the  axis  corresponding  to  Xi(y)  from  point  ra* . 

We  also  must  search  for  the  best  value  of  h  graphically.  In  order  to  help  find  the  best  h,  we  plot 
Xi(y)  versus  both  X2(y)  and  x^OO,  where  the  latter  is  defined  as 

X'2(y)  =  vi/Ntot(l-y)  ],  (6.31) 

and  look  for  a  close  fit.  Then,  we  define  H'  as  the  convex  hull  of  the  plot  of  x'2  versus  Xi ,  and  m'2 
as  the  moment  of  x'2(Ii)  •  Note  that  the  intersection  of  H  and  H'  contains  moments  corresponding  to 
symmetric  probability  distribution  functions  in  P  ( 2 ).  Moreover,  the  convex  set  H  for  h  is  the  same 
as  the  set  H'  for-h.  Therefore,  we  only  need  to  consider  positive  values  of  h.  (  Note  that  m2  =  m'2  .) 

For  plotting  x2(y)  and  xVy)  versus  Xi(y) .  we  need  to  evaluate  4^  for  values  of  y  in  the  range 
E.  For  the  e-mixture  example,  we  use  (6.12)  and  draw  the  curve  for  various  distribution  functions  in 
this  model.  Figures  6.2  through  6.7  depict  the  sets  H  and  H’  for  these  distributions  and  for  N  -  31 


and  63  with  the  same  DS/SSMA  system  configurations  considered  in  previous  examples.  Notice  that 
in  these  plots  the  convexity  of  these  sets  is  not  apparent  since  the  vertical  axes  are  in  logarithmic 
scale  and  the  horizontal  axes  are  in  linear  scale.  We  apply  the  moment-space  bound  to  two  and  six- 
user  DS/SSMA  systems  in  impulsive  channels.  The  results  are  summarized  in  Tables  6.2  and  6.3.  As 
expected,  for  a  fixed  number  of  users  K,  as  N  increases,  D  decreases  and  the  bounds  tend  to  become 
tighter.  On  the  other  hand,  with  fixed  N,  as  the  number  of  users  increases,  the  bounds  become  looser. 
Similar  to  what  happened  for  the  Taylor-series  approximation,  the  moment-space  bound  results  in 


better  performance  approximation  for  more  impulsive  channels.  A  possible  explanation  for  this  is 
that  the  exponential  kernel  in  (6.25)  better  approximates  X2(y)  given  in  (6.23)  for  impulsive  channels 
than  it  does  for  the  Gaussian  case.  (  Note  from  Table  6.1  that  the  Taylor  series  approximation  lies 


III 

\m 

vv 

It'**-  I 


within  the  moment-space  bounds  in  all  cases  common  to  both  sets  of  results. ) 


6.4.  Summary 


In  this  chapter  we  have  used  two  techniques  to  estimate  the  average  error  probability  of  the 
linear  correlation  receiver  in  non-Gaussian  impulsive  channels.  These  methods  were  introduced  as  an 


it 

m 


! 


alternative  to  the  exact  computation  proposed  in  Chapter  4  when  N  or  K  is  large.  Particularly  for 
impulsive  channels,  both  techniques  appear  to  estimate  the  performance  closely.  The  principal 
difficulty  that  arises  in  applying  these  methods  is  that  the  probability  distribution  function  of  the  sum 
of  N  non-Gaussian  random  variables  has  to  be  evaluated  for  some  points  near  the  signal-to-noise 
ratio.  However,  this  problem  is  essentially  that  of  single -user  analysis,  and  for  the  e-mixture  case  it  is 


trivial. 
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TABLE  6.2.  MOMENT -SPACE  UPPER  AND  LOWER  ERROR  BOUNDS 
FOR  LINEAR  CORRELATION  RECEIVERS  IN  THE  BINARY  PSK 
DS/SSMA  SYSTEM  ;  GAUSSIAN  AND  e-MIXTURE  CHANNELS,  SNR-8.0 
dB,  N-31. 


DISTRIBUTION 

K 

P 

— - -« - — 

K 

2 

3.93  x  1(T* 

4.16  xlO-4 

4.64  x  ir4 

GAUSSIAN 

6 

1.44  x  10-4 

3.48  x  ir3 

8  .58  x  1(T3 

e-MIXTURE 

2 

3.41  x  10"3 

3.46  xl(T3 

3.47  x  10"3 

e-0.01,  f-100 

6 

4.10  xir3 

8.80  xlO"3 

2 

1.31  x  ir3 

1.31  x  1(T3 

1.35  xir3 

e-O.l.y-lOO 

6 

1.51  x  1<T3 

8.80  x  ir3 
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TABLE  6.3.  MOMENT- SPACE  UPPER  AND  LOWER  ERROR  BOUNDS 
FOR  LINEAR  CORRELATION  RECEIVERS  IN  THE  BINARY  PSK 
DS/SSMA  SYSTEM  ;  GAUSSIAN  AND  e-MIXTURE  CHANNELS,  SNR-8.0 
dB,  N=63. 


DISTRIBUTION 


GAUSSIAN 


e-MIXTURE 
e-0.01,  -^-100 

e-0.1,  t2— 100 


2 

2.96  x  KT* 

3.06  x  KT* 

3.24  xHT* 

6 

3.91  x  KT* 

1.33  x  1(T3 

2 

1.72  xl(T3 

1.76  x  1(T3 

1.77  x  1(T3 

6 

2.13  x  kt3 

3.39  x  1(T3 

2 

7.19  xl(T3 

7.19  x  1(T3 

722  x  1(T3 

6 

8.97  x  KT4 

2.48  x  1(T3 

m 
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CHAPTER  7 


HARD-LIMITING  CORRELATION  RECEIVERS  IN  IMPULSIVE  CHANNELS: 


APPROXIMATIONS 


One  conclusion  from  the  hard-limiting  correlation  receiver  analysis  in  Chapter  5  is  that  comput¬ 


ing  the  exact  expression  for  Pe  can  be  very  time  consuming  for  large  N.  For  this  reason,  when  signa¬ 


ture  sequences  of  length  more  than  31  are  used,  bounds  and  approximations  are  needed  to  replace  the 


exact  computation  of  Pe.  In  this  chapter  we  propose  an  approximation,  an  upper  bound,  and  a  lower 


bound  for  the  multi-user  error  probability  of  the  hard-limiting  correlation  receiver  in  impulsive  noise. 


The  lower  bound  uses  the  idea  of  truncating  a  finite  sum  of  positive  numbers  in  the  expression  for  the 


error  probability  (5.18);  the  approximation  is  based  on  estimating  the  probability  distribution  function 


of  the  test  statistic  YN;  and  the  upper  bound  uses  a  Chemoff  bounding  technique.  Also  in  this 


chapter  we  analyze  the  hard-limiting  correlation  receiver’s  performance  asymptotically  in  N. 


7.1.  A  Truncated  Series  Lower  Bound  on  Pe 


Our  first  bound  on  the  error  probability  of  the  hard-limiting  correlation  receiver  uses  the 


property  that  Pr  YN  =  2  m- 1  |  b(JI)= — 1,  pi  decreases  sharply  as  m  increases  to  (N  +  l)/2.  The 


event  {  Y^  =  m  }  has  very  small  probability  for  large  values  of  m  when  conditioned  on  b<^=— 1, 


because  YN  taking  some  positive  value  m  means  the  output  of  the  hard-limiter  was  +1  exactly 


(N  +  m)/2  times,  an  event  which  becomes  much  less  likely  as  m  increases.  Therefore,  in  computing 


Pe  via  (  see  also  (5.18) ) 


(N+O/2  i  2  r 

P.CD*  I  TrZ  /  Pr  YN  =  2m-l  |  bjn 

m«l  b  0  L 


=“l.t,  b,  <i>  d<t», 
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we  can  obtain  a  reasonably  tight  lower  bound  by  adding  only  the  first  few  terms  of  the  sum  of  odd 
m’s  from  1  to  N.  The  accuracy  of  this  approximation  is  easily  controlled  by  including  more  terms  in 
the  computation,  noting  the  tradeoff  with  the  cost.  As  an  example,  we  compute  a  lower  bound  by 
including  only  |^log2(N+l)j-l  terms  in  the  summation  in  (7.1)  instead  of  (N+l)/2.  For  the  e- 

mixture  and  Laplacian  examples,  we  examine  the  tightness  of  this  lower  bound  of  Pe(T)  for  the 
hard-limiting  correlation  receiver.  Typical  results  showing  the  accuracy  of  this  approximation  are 
depicted  in  Figure  7.1  and  Tables  7.1  and  7.2.  (  Also  shown  are  approximate  values  to  be  discussed 
below.  )  These  results  are  obtained  for  the  same  DS/SSMA  system  considered  in  examples  in 
Chapter  5.  Examples  are  carried  out  with  lengths  of  signature  sequences  equal  to  31  and  63.  Tables 
7.1  and  7.2  correspond  to  SNR  =  8  dB  and  SNR  =  4  dB,  respectively.  As  we  see  from  these  tables 
and  the  figure,  the  bound  is  generally  tight.  Moreover,  for  more  impulsive  examples,  which  are  the 
most  interesting  ones,  the  lower  bound  is  extremely  tight. 


7.2.  Binomial  Approximation 

The  above  bound  is  still  computationally  expensive,  although  much  less  so  than  exact  computa¬ 
tion.  Thus,  simpler  estimates  of  Pc  are  of  interest.  One  such  estimate  can  be  based  on  approximating 
the  distribution  function  of  the  test  statistic  YN ,  which  is  given  as 


M-i  N-i 

YN=  Z  Zj(1)=  £  sgn 


j=o 


j=0 


(7.2) 


The  error  probability  for  the  hard-limiting  correlation  receiver  can  be  written  in  a  form  slightly  dif¬ 
ferent  from  (7.1),  namely, 


*4 


i 

r*mT 


& 


*' 

$ 

'V, 


3 

£ 


85 


TABLE  7.1.  TRUNCATED  SERIES  LOWER  BOUND  AND  BINOMIAL 
APPROXIMATION  FOR  THE  ERROR  PROBABILITY  OF  HARD-LIMITING 
CORRELATION  RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ; 
GAUSSIAN  AND  IMPULSIVE  CHANNELS,  SNR-8.0  dB,  K-2  AND  TYPI¬ 
CAL  TIME  DELAYS. 


DISTRIBUTION 

N 

IHK3HI 

P.= 

P. 

31 

4.51  x  ir3 

6.92  xlOr3 

4.63  x  ir3 

GAUSSIAN 

63 

2.45  x  1CT3 

2.93  x  1<T3 

2.51  x  10r3 

31 

8.55  xHT* 

1.77  x  ir3 

8.67  x  10^ 

LAFLACIAN 

63 

1.01  x  104 

1.63  x  10T4 

1.02  x  10^ 

e-MIXTURE 

31 

3.87  x  10-4 

1.16  xKT3 

3.90x10“* 

£-0.01, -y2- 100 

63 

5.73  x  10rs 

1.05  x  KT* 

5.77  x  101*5 

31 

2.58  x  10"* 

4 .54  x  1CT5 

2.59  x  lO-* 

6—0.1,  y2— 100 

63 

2.06  x  10"10 

1.68  xl(T8 

2.06  x  10'10 

TABLE  7.2.  TRUNCATED  SERIES  LOWER  BOUND  AND  BINOMIAL 
APPROXIMATION  FOR  THE  ERROR  PROBABILITY  OF  HARD-LIMITING 
CORRELATION  RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ; 
GAUSSIAN  AND  IMPULSIVE  CHANNELS,  SNR-4.0  dB,  K-2  AND  TYPI¬ 
CAL  TIME  DELAYS. 


DISTRIBUTION 

N 

BKSHI 

P. 

31 

3.89  xir2 

4.52  x  1<T2 

4.14x10 r2 

GAUSSIAN 

63 

3.45  x  ir2 

3.81  x  ir2 

3.68  xir2 

31 

9.77  x  1(T3 

1.26  xir2 

1.01  x  ir2 

LAPLACIAN 

63 

4.70  x  1(T3 

5.37  x  1(T3 

4.84  x  10*3 

e-MDCTURE 

31 

9.39  xl(T3 

1.27  x  ir2 

9.72  xir3 

e-0.01,  100 

63 

6.33  x  10"3 

7.21  x  10'3 

6.54  xir3 

31 

6.63  x  1(T5 

3.42  x  IQT* 

6.66  x  10'5 

e-0.1,  7**100 

63 

7.60  x  10-7 

3.78  x  10-6 

7.61  x  10*7 
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Pe(S)  =  Eb,«  '  Pr  I  bo=+l  ,  p/1J  j  j*  . 


Note  that,  conditioned  on  the  multiple-access  parameters,  the  test  statistic  (Ym+N)/2  is  multinomially 
distributed.  Here  we  approximate  this  distribution  by  a  binomial  distribution  for  large  N  because  the 


approximation  is  asymptotically  exact  as  N  grows  without  bound,  and  the  binomial  distribution  has 


many  attractive  analytical  properties.  Thus,  we  approximate  (7.3)  as 


Pe(l)  =  Eb.J  [Njq/(l-qi)N-j  , 


where  qx  is  the  "success"  probability  associated  with  the  approximating  binomial  distribution.  To 


choose  qi  we  equate  the  mean  of  the  binomial  approximation, 


E^YnI^+I.p,  f=N[  l-2qj], 


with  the  actual  mean  of  the  test  statistic. 


e|yn  |btf>-+l,Pl  -=N-2N£  Prj^l+l/”] 


For  the  e-mixture  example  (  which  includes  Gaussian  noise  as  a  special  case  )  qj  is  given  as 


i  n-i  r  fi+ij(1)]  f  i+ij(i) 

<ii  =  TT  2  0-*)Q  -d-  -d- 
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where 

(1-e+ey2)  2Eb(1) 


For  the  Laplacian  example,  qi  is  computed  from  die  following 


expression 


(7.8) 


.  ,A  NoN 
where  <r  =  — — , 
2E^> 


and  we  have  assumed  that  |  Ij(1)  |  S 1,  which  holds  for  example  when  K  -  2  or 


when  e^i  -  1  for  all  k.  When  K-2,  (7.7)  and  (7.8)  can  be  simplified  by  noting  that  |  Ij(1)  |  takes  only 
two  forms  as  in  (5.15).  Substituting  (7.7)  or  (7.8)  into  (7.4)  yields  the  binomial  approximation  for  the 
error  probability  Pe(T )  for  the  hard-limiting  correlation  receiver. 


Figure  7.1  and  Tables  7.1  and  7.2  contain  some  numerical  results  that  compare  this  binomial 
approximation  to  the  above  bound  and  to  the  exact  values  of  Pe(T)  for  different  values  of  N.  (  A 
DS/SSMA  system  identical  to  the  one  used  in  previous  chapters  is  considered. )  The  approximation 
is  evaluated  for  various  parameters  of  e-mixture  noise  and  also  for  Laplacian  noise.  Considering  the 
simplicity  of  the  binomial  approximation,  the  results  in  Figure  7.1  and  Tables  7.1  and  7.2,  particu¬ 
larly  for  large  N,  indicate  that  (7.4)  is  a  reasonably  close  estimate  of  Pe(T).  In  Section  7.4.  we  will 
see  that  if  signature  sequences  are  chosen  appropriately,  the  binomial  approximation  is  asymptoti¬ 
cally  exact 


7.3.  Chernoff  Bound 

In  this  section,  we  examine  the  tightness  of  the  Chernoff  bound  when  applied  to  the  hard- 
limiting  correlation  receiver  problem.  Saltzberg  in  [30]  and  Lugannani  in  [15]  used  this  method  for 
intersymbol  interference  error  bounding,  which  is  similar  to  the  problem  of  interest  here.  They  found 
the  Chernoff  bound  to  be  superior  to  existing  upper  bounds  in  intersymbol  interference  studies.  How- 


ever,  our  investigation  (  not  reported  in  this  thesis  )  showed  that  applying  the  Chemoff  bound  to  the 


multi-user  communication  system  in  impulsive  noise  with  the  linear  correlation  receiver  results  in  a 


loose  bound  for  error  probability.  For  the  linear  case  the  moment-space  bounds  are  better.  However, 


the  moment-space  technique  cannot  be  applied  directly  to  the  hard-limiting  receiver,  and  so  despite 


the  possibility  of  achieving  a  loose  bound,  we  pursue  the  Chemoff  technique  for  the  hard-limiting 


correlation  receiver. 


To  understand  the  application  of  the  Chemoff  bound,  we  first  write 


Pe(l)  =  E^.b|pr  [ynS1  Ib^W— 1,Pi]|. 


Then,  we  define  a  new  zero  mean  random  variable  as 


f\j  =  Zj(I)  —  E  {  Zj(1)  |  =-l ,  pi }  , 


(7.10) 


where  Zj(1)  is  given  in  (7.2).  Now  we  can  rewrite  (7.9)  as 


Pe(I)  =  E0  jPr  "SVI-  ZMb^—l.pi  [. 

L  i-o  i-o  JJ 


(7.11) 


where 


Uj  =  E  |  ±/l)  |btf>— l,p,  1=24^  -I 


(7.12) 


In  (7.12)  'F(a)  is  defined  as  Pr[rj>a  ],  where  r\  is  a  typical  normalized  noise  sample  with  unit  vari¬ 


ance.  Now  if  we  define  p.s  as  £  )ij  and  then  we  can  write 

j-o 


r*;i 


90 


Pe(l)  =  E*,b 


N+l 

2 


Z  Pr  [f|,  =  2m-l  |  bil)=-l,pIl 

B-l  L  J 


(7.13) 


an-i 

where  f),  =  X  Aj-  Now,  for  each  t  >  0 ,  we  introduce  a  discrete  random  variable  with  mass  func- 

j«0 


non 


Pr[Ct  =  a]  = 


euPr[f>t==a] 


,  t>0 


(7.14) 


where  is  the  moment  generating  function  of  the  random  variable  f|s.  Replacing  f),  with  £t>  the 
expression  inside  the  curly  bracket  in  (7.13)  can  be  bounded  from  above  as 


N+I 

2 


Z  Pr  [f)s  =  2m-l-  JaM^ft)  X  Pr  k  =  2m-l-m] 
n»l  L  m»  l  ^  J 


»t((i,+l-2m) 


N+I 


^  exp  [m^(t)  -  t(l— jo.*)  )  X  Pr  Ci  =  2m-l-n,  ,  (7-15) 

m-1 


where  m^ft)  =  In  [  (t)  ].  Finally,  for  any  value  of  t  on  the  positive  real  line  we  have 


(7.16) 
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V 


To  get  the  best  bound  we  minimize  the  right-hand  side  of  (7.16)  with  respect  to  t  In  particular, 
we  are  looking  for  to  such  that 


to  =  arg 


(7.17) 


Differentiating  the  argument  in  (7. 17)  and  setting  the  result  equal  to  zero  results  in 


m'fli(t)=  1-ix,  . 


(7.18) 


In  order  to  see  whether  a  solution  to  (7.17)  exists  and  it  satisfies  (7.18),  we  write  the  expression  for 
m^t  (t)  for  the  hard-limiting  correlation  receiver  as 


N-l 

*V1)*  £  In 

j»0 


(7.19) 


Using  (7.19),  it  is  easy  to  show  that  to  satisfying  (7.18)  is  obtained  by  solving 


n-i  i-nj-(l+Mj)e2t» 
i-o  1-Pj  +  (1+Pj)e2t0 


=  0  ,  to>0  . 


(7.20) 


From  (7.19)  it  follows  that  a  solution  to  (7.20)  exists  and  satisfies  (7.17).  Finally,  the  Chemoff 
bound  on  the  error  probability  of  the  hard-limiting  correlation  receiver  is  given  as 


Pe(l)SEb,« 


[to Us  -to  +  (to) 


(7.21) 


92 


where  to  is  obtained  by  solving  (7.20).  We  examine  the  tightness  of  the  Chemoff  upper  bound  on 
Pe(t)  for  the  e-mixture  and  Laplacian  examples.  Typical  results  accessing  the  accuracy  of  this 
approximation  are  depicted  in  Tables  7.3  and  7.4.  These  results  are  obtained  for  the  same  DS/SSMA 
system  considered  in  the  examples  in  Chapter  5.  These  examples  are  carried  out  with  lengths  of  sig¬ 
nature  sequences  equal  to  31,  63,  127.  Tables  7.3  and  7.4  correspond  to  SNR  -  8  and  SNR  -  4, 
respectively.  As  expected,  this  bound  is  not  particularly  tight;  however,  its  simplicity  and  the  lack  of 
alternative  bounds  makes  it  potentially  useful. 


7.4.  Asymptotic  Analysis 

In  Section  7.2  an  approximation  was  obtained  for  the  bit-error  probability  of  hard-limiting 
correlation  receivers  in  impulsive  noise  that  was  good  for  large  N.  Also  in  Chapter  5,  more  significant 


improvement  in  performance  of  the  hard-limiting  correlator  over  the  linear  correlator  was  experi¬ 
enced  as  we  used  longer  signature  sequence  lengths.  Motivated  by  these  findings,  here  we  investigate 
the  behavior  of  hard-limiting  DS/SSMA  correlation  receivers  when  the  length  of  the  signature 
sequences  used  per  data  bit  becomes  infinitely  large. 

Recall  from  the  asymptotic  analysis  in  Section  6.2  that  conditioned  on  the  multiple-access 
parameters  pb  the  test  statistic  is  a  sum  of  N  independent  random  variables.  Therefore,  assuming  that 


E{YNlbjlWl,p,} 

/  Var(  Yn  !  b^Wl.p,) 


exists  almost  surely , 


(7.22) 


the  as-erage  bit-error  probability  can  be  written  asymptotically  as 

.  _  LfEvll 


Q  o7 1  f ■ 


(7.23) 
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TABLE  7.4.  CHERNOFF  UPPER  BOUND  FOR  THE  ERROR  PROBABILITY 
OF  HARD-LIMITING  CORRELATION  RECEIVERS  IN  THE  BINARY  PSK 
DS/SSMA  SYSTEM  ;  GAUSSIAN  AND  IMPULSIVE  CHANNELS,  SNR-4.0 
dB,  K-2  AND  TYPICAL  TIME  DELAYS. 
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where  Ey  =  lim  E  ■{ 

N-*+~  1 


-7=  |  b^’Wl  ,p|l  and  (Oy)2=  lim  Var 


.  Note  that 


•vN  ’  '  '  '  I  '  '  N- 

here  is  given  in  (7.2).  For  a  given  set  of  multiple-access  parameters  the  normalized  mean  can  be 
written  as 


1  N-l 

Ey  =  lim  4=  X  ^ 
N-*+-  VN  jto  I 


1-2T 


(7.24) 


where  'F(a)  =  Pr  [  r(  >  a  ]  is  the  complementary  distribution  function  of  a  typical  noise  sample  rj,  nor¬ 
malized  to  have  zero  mean  and  unit  variance.  For  the  e-mixture  example,  the  distribution  function  in 
(7.24)  is  given  as 


'F(a)  =  (l-e)Q  [aVl-e+eT2]  +eQ  a^e+^p 


(7.25) 


Similariy  for  the  Laplacian  example  we  have 


*F(a)  =  ~  exp  {  -V2a  }  ,  fora>0 


(7.26) 


Referring  to  (7.24),  the  probability  distribution  function  can  be  expanded  in  Taylor  series  about 


a 


VN 


as 


'F 


Jrti+i. /»]  . 


VN '  “  JL 

n*  1  ^ 

N2n! 


VN’ 


(7.27) 


where  vF(nl  is  the  riJ*  derivative  of  'V.  Now  if  H'C.)  is  analytic  and  we  have  the  following  asymptotic 


-.7  V 

•>3 


•  *  *  1 

.  .  . . . % ^ 


(l  '1  ^ 

vi 


I 
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‘>•1 

*y.»1 


.’M 


52a 

m 

mmH 


♦S« 


property: 


N-l  ~  an(Ii(1))n  rt 
j*u  n»i  2  , 

N  *  n! 


(7.28) 


then  we  can  write 


Ev=NUm_VN  1-24-1^)  -2a^(0)  , 


(7.29) 


where  is  the  probability  density  function  of  rj .  From  (7.2)  and  (7.29),  it  follows  that  ( cry  )2  is 


equal  to  1.  This  yields  a  simple  expression  for  the  asymptotic  error  probability  of  the  hard-limiting 


correlation  receiver  as 


N!™„pe=Q  12a%(0) 


(7.30) 


if  (7.28)  is  satisfied.  Note  that  Q(2  af^(0) )  is  the  asymptotic  single-user  performance  of  the  hard- 


limiting  correlator  (  see  Chapter  3  ).  Thus,  (7.30)  implies  that  single-user  asymptotic  performance  is 


achieved  by  the  hard-limiting  correlation  receiver  if  the  condition  (7.28)  is  satisfied. 


Regarding  condition  (7.28),  note  that  the  order  of  summations  can  be  interchanged  if 


n-i  -  an(Ij(1))n 


2  I  -1ST— <+~  <»•»•> 

N  — >  +  «  •  q  n* i  — 

N  2  n! 


(7.31) 


Assuming  that  |l.(l)|<K  and  that  there  exists  ueR  such  that  I  4^nV  —£= )  <u  for  all  n  and  large 

\N 


enough  N  then  condition  (7.31)  is  satisfied  if  the  number  of  users  is  finite.  However,  (7.31)  can  hold 
for  infinitely  large  K,  for  instance,  it  is  easy  to  see  that  (7.31)  will  be  satisfied  if 


N-i  |l|" 

£  -1tr<+°’  <a  s->  v”  - 
i-o  .,-r 


By  induction  we  can  show  that  the  conditions 


N-l  |IP|2 

lim  I|(b,T,<J»)  =0  and  lim  V  — rr-  =  0  (a.s.)  , 

N  — *+«•  —  —  N  — »+•• 


(7.32) 


(7.33) 


are  sufficient  to  meet  (7.28)  and  (7.31)  if  K£N1/4.  Note  that  the  first  part  of  (7.33)  is  identical  to  the 
condition  for  the  linear  correlator  (6.21).  The  second  part  is  satisfied  by  imposing  slighdy  more  res¬ 
trictive  upper  bound  on  the  number  of  users.  Thus,  the  ability  to  achieve  a  single-user  performance 
asymptotically  in  N  is  possible  within  the  above  regularity  conditions,  and  so  Table  3. 1  and  Figures 
3.8  and  3.9  contain  the  asymptotic  average  bit-error  probability  for  the  hard-limiting  correlation 
receiver  within  these  conditions. 
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CHAPTER  8 

NONLINEAR  CORRELATION  RECEIVERS  IN  IMPULSIVE  CHANNELS: 

ASYMPTOTIC  ANALYSIS 

In  Chapters  6  and  7,  we  investigated  the  asymptotic  performance  of  linear  and  hard-limiting 
DS/SSMA  correlation  receivers  .as  the  lengths  of  the  signature  sequences  grew  without  bound.  In 
turn,  asymptotically  ideal  DS/SSMA  signals  were  proposed.  To  gain  further  insight  into  die  non¬ 
linear  problem  and  explore  the  fundamental  limitations  of  these  systems,  we  carried  out  an  asymp¬ 
totic  analysis  for  general  nonlinear  DS/SSMA  correlation  receivers  depicted  in  Figure  8.1,  where  g  is 
a  general  memoryless,  nonlinear  element.  Then,  we  searched  for  necessary  asymptotic  properties  for 
the  spreading  sequences  that  assure  asymptotic  single-user  performance  in  a  multi-user  environment. 

Recall  that  in  the  asymptotic  analysis  the  energy  per  data  bit  for  a  user,  the  length  of  the  data 
period  T,  the  noise  power,  and  in  turn  the  signal-to-noise  ratio  are  kept  constant  However,  the  length 
of  the  signature  sequences  used  per  data  bit  is  growing  without  bound,  and  in  turn  the  sampling  (  and 
chip  )  period  Tc  is  decreasing  to  zero.  Also  note  that  we  rescaled  the  data  by  VN  to  keep  the  correct 
scale  (  this  prevents  us  from  having  to  allow  the  nonlinearity  to  change  with  N  ),  and  therefore  the 
test  statistic  is  written  as 


N-l 

Yn=  I  g 

j  =  0 


(8.1) 


N0T 

where  T]j ’s  are  samples  of  noise  with  zero  mean  and  variance  —  and  g  is  the  nonlinear  element.  In 

(8.1)  the  I'jW’s  denote  the  contribution  of  multiple-access  interference  in  the  samples  of  the  received 
signal  and  are  defined  as 


j=0.1 . N-l, 


(8.2) 


where  Ij(1)  is  given  in  (5.12).  Note  that  l'[l)  relates  to  Ij(b,T,<t>)  defined  earlier  in  (4.5)  through 


Ii(b,x,«)  = 


(8.3) 


Referring  to  (8. 1),  conditioned  on  the  multiple-access  parameters  pi,  YN  is  seen  to  be  the  sum  of 
N  independent  random  variables.  To  apply  the  central  limit  theorem  to  the  sum  (8.1),  we  first  assume 
that 


E  {  Yn  1  bj^=+l ,  pi } 

N  -*  +  «*»  VVar(YN|bi»»+ltPi) 


exists  almost  surely , 


(8.4) 


then,  via  the  bounded  convergence  theorem,  the  average  bit-error  probability  for  the  nonlinear 
DS/SSMA  receiver  can  be  written  asymptotically  as 


1 


8 


8 

1 

1 


1 


Pe  —  Ej.o.b 


(8.5) 


where  Ev  =  lim  E  I  — ==■  |  b(P=+l , p,  >  and  (ay)2  =  lim  Var 
N— ♦  +  «•  !  vN  I  N— >+•» 


For  a 


given  set  of  multiple-access  parameters  the  normalized  conditional  mean  can  be  written  as 
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E?  =  N!£j-  W  [s(x  +  0N[l  +  l/1)])-g(x-0N[l+l/1)])]fli(x)dx,  (8.6) 


where  0n  =  V^T.  Similarly,  the  variance  of  the  normalized  test  statistic  has  the  form 


( oy )2  =  Nlhn_  V  f  [g2  ( x  +  0N  [1  + !/» ] )  +  g2  ( x-0N  [1  + 1/1)  1 )  ]  f^(x)  dx 


l  N-l 

■m55-n  £ 


J  [g(x+eN[i +ij(,)])-g(x-eN[i+ij(1)]) j  f^(x>dx 


(8.7) 


In  (8.6)  and  (8.7)  if  conditions  of  the  dominated  convergence  theorem  are  met  by  the  integrands, 
then  the  order  of  the  infinite  sums  and  the  integrals  can  be  interchanged  (see  [6]  page  257).  We  first 
concentrate  on  (8.6);  if  g(-)  is  analytic  and  has  a  bounded  second  derivative  then  we  expand  the 
integrands  in  (8.6)  as 


0M2  I  I(1)  I2 

g(x  +  0N[l  +  l/1J  ])  =  g(x  +  0N)  +  0N  g'  (  X  +  0n  )  H - — — ~ - —  g"  (  £  ), 


(8.8) 


for  some  ^e(x  +  0N,x  +  0N  [  1+I/I)] ).  Substituting  (8.8)  into  (8.6),  we  realize  that  if  we  have  the  fol¬ 
lowing  asymptotic  properties 


sup(g")  lim  —  £  jlj(1)j2  =  0  (  a.s. ), 


3 

"2  i*° 


(8.9) 
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N-ll/O  +!*  r  1 

J  [g,(x  +  eN)  +  g'(x-eN)jf1i(x)dx  =  0  (a.s.). 


(8.10) 


j«0  ^  0 


A  ^  /  2Eb(1)  nA^n^T 

with  a=  — — ,  and  p=\f  — - — ,  then  we  can  simplify  (8.6)  as 


+  «*  - 

Ey=  J  Jim  VN  g(x  +  0N)-g(x-0N)  f^(x)dx, 
n  N-*+»  J 


(8.11) 


where  f^(.)  is  the  probability  density  function  of  a  typical  zero  mean  noise  sample  which  has  variance 
P2.  Note  that  (8. 1 1)  can  be  simplified  further  to 


Ey=  Japg,(x)f1i(x)dx, 


(8.12) 


Similarly  for  (ay)2  in  (8.7),  since  g2  is  analytic  and  assuming  that  g2  has  bounded  second 
derivative,  then  the  integrands  can  be  expanded  in  a  Taylor  series  about  a  point  near  x,  i.e., 


On.2  1 1(1>  l2 

g2(x  +  0N[l  +  l/,>])  =  g2(x  +  0N)  +  0Nl/1>  (g2)'(x  +  0N)  +  -^— y — —  (g2)"^),  <8-13) 

for  some  i;e  ^x+0N,x  +  0N[  1+Ij(1)]j.  In  (8.13)  (g2)"  is  the  2nd  derivative  of  g2.  Substituting 
(8.13)  into  (8.7),  we  realize  that  if  we  have  the  following  asymptotic  properties 


-^-sup((g2)")  Jim  SI|I/,)I2  =  0  (a.s.), 

2  N->+«  Nz  j_0 


(8.14) 


i 

1 
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aPN!^«~T "jVj  [(g2),(x  +  0N)-(g2),(x-0N)]f1i(x)dx  =  O  ( a . s . ).  (8.i5) 

N2  J” 


and  also 


g(a  +  eN[l  +  l/l)])-g(x- 


(x)dx 


(a.s.)  , 


(8.16) 


for  some  d  >  0  and  all  j  »  0, 1 . N  - 1,  then  we  can  simplify  (8.9)  as 


( )2  =  2  J  g2  (x)  (x)  dx. 

o 


(8.17) 


Finally,  (8.12)  and  (8.17)  yield  the  following  expression  for  the  argument  of  the  Q  function  in 

(8.5): 


E V 
oy 


2<*P  J  g'  (x)  fvj  (x)  dx 


2jg2(x)fTj(x)  dx 


(8.18) 


if  the  conditions  (8.9),  (8.10),  (8.14),  (8.15),  and  (8.16)  are  satisfied.  By  substituting  (8.18)  into  the 
expression  for  the  asymptotic  error  probability  given  in  (8.5),  we  recognize  the  possibility  of  single- 
user  performance  as 


lim  Pe  =  Q 

N  — *  +  «> 


(8.19) 
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where  v  is  the  detection  efficacy  defined  in  (3.11),  which  for  the  nonlinear  correlation  receiver  is 
easily  derived  as 


v  = 


}g'(x)f^(x)dx 


Jg^xKnOO  d* 


(8.20) 


We  have  just  shown  that  asymptotically  in  N  and  within  regularity  on  the  nonlinearity  the  non¬ 
linear  correlation  receiver  can  perform  with  single-user  error  probability  in  a  multi-user  environment 
if  the  conditions  (8.9),  (8.10),  (8.14),  (8.15),  and  (8.16)  a  j  satisfied.  Since  we  assumed  that  both  g 
and  g2  have  bounded  second  derivative,  the  above  conditions  are  satisfied  if 


lim  ^(b.x,^)  =0 

N  -*+~  “  — 


and 


lim 

N-*+« 


N-J 


£ 

j-0 


(a.s.). 


(8.21) 


if  the  number  of  users  does  not  grow  faster  than  NI/4.  Note  that  the  first  part  of  (8.21)  is  identical  to 
the  condition  for  the  linear  correlator  (6.21).  The  second  part  is  satisfied  by  imposing  slightly  more 
restrictive  upper  bounds  on  the  number  of  users,  and  is  the  same  as  the  analogous  condition  derived 
for  the  hard-limiting  correlator. 

To  study  further  the  condition  (8.21),  we  first  recall  from  (6.2)  that  the  multiple-access  interfer¬ 
ence  can  be  written  as 


I 


K  -  Wkl 

Ii(b,T,<t>)=  £  cos<t>k  — 7- 

lc-2  ^ 


(8.22)  ■ 

I 


where 
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WU  =  bi?  [ dkCM  ( mk  -  N  +  1)  +  (l-dk)  Cu  ( mk  -  N )  ] 


+  b^  [dkC^i  (mk+l)  +  (1-dfc)  C^i  (m^)]. 


(8.23) 


A  Ntk 

In  (8.23)  dk  =  - 


NX, 

A 

Nxk 

T 

,mk  = 

T 

L 

,  and  C^j  is  the  discrete  aperiodic  cross-correlation 


function  defined  in  (4.10).  For  a  finite  number  of  users  and  from  (8.22)  and  (8.23)  it  follows  that 
(8.21)  is  satisfied  if  we  have  the  following: 


n3S-  n  "2*  [c^  <”)]=»• 


(8.24) 


for  all  k,i  such  that  k*i.  This  is  equivalent  to  the  conditions 

lim  max  ■—  |  O^m)  |  =0  and  lim  max  —  I  0kj(m)  1=0, 

N  —>+■»  m  IN  ra  N 


(8.25) 


for  all  k  and  i  such  that  k*i.  As  mentioned  earlier,  the  existence  of  binary  sequences  satisfying  these 
conditions  follows  by  applying  (8.24)  to  infinite  sequences  proposed  in  [34],  For  these  sequences  the 
convergence  of  the  crosscorrelations  to  zero  is  very  slow.  However,  for  the  infinite-user  case  the  con- 

l 


dition  (8.22)  is  still  satisfied  as  long  as  the  number  users  do  not  grow  faster  than  N4 . 

An  important  issue  to  discuss  here  is  the  interpretation  of  the  desired  condition  (8.25)  in  terms 
of  bandwidth  requirements.  First  we  write  the  discrete  Fourier  transform  of  the  signature  sequence 


a(k)  as 


A  1  N"> 


Am=^r  Saf'WrT  m=0,l . N- 1, 

w  i-o 


(8.26) 


ot 


V 

v 


where  WN  =  exp  {  -i-^- }.  Let  ©iy(0  denote  the  DFT  of  the  periodic  crossconelation  function  0y(-). 
It  is  well  known  that 


0y(m)  =  NA.m«Ara®  m=0,l . N-l. 


(8-27) 


From  the  discussion  in  [34],  a  large  set  of  infinitely  long  sequences  can  be  generated  such  that  for 
these  sequences  (8.25)  approaches  zero  as  fast  as  N-11  for  some  0<u<-—.  For  the  sequences  in  the  set 

Sn.  ©kj(ra)  grows  with  N,  but  as  slow  as  N1_u  with  0<u<-^-  for  all  m.  From  (8.27)  it  follows  that 

A-m*  Am(,)  approaches  zero  as  fast  as  N-u  for  all  m  and  all  pairs  of  sequences  in  Sn-  In  summary, 
(8.25)  implies  that 


lim  A.,,,®  Am(l)  *  0  for  all  meZ. 

N-»+— 


(8.28) 


This  indicates  that  the  signals  are  occupying  different  frequency  bands  which  causes  the  infinite-user 
communication  system  with  asymptotically  ideal  sequences  to  use  infinitely  large  bandwidth. 

To  illustrate  the  frequency  domain  characteristics  of  these  sequences  consider  the  following 
example.  Suppose  the  signature  sequence  a(1)  with  a/l)=  +  1  for  all  0<j<N  - 1  is  a  member  of  the  set 
of  asymptotically  ideal  sequences.  Then  the  DFT  of  a(1)  is 


A  ri)_ 

rVm 


1  if  m=0 

0  if  m=l,  2 . N-l 


(8.29) 


Now  as  N  grows  without  bound,  the  condition  (8.28)  impiies  that  all  other  signature  sequences  in  SN 
must  have 
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Hypothetically,  the  set  of  K  sequences  can  be  built  by  forcing  the  DFT  of  a(k)  to  have  a  nonzero  value 
only  at  one  point,  A j[k\  *  0.  Thus  (8.28)  is  satisfied;  however  the  spectrum  is  completely  occupied  by 
the  infinitely  many  users.  This  is  not  surprising,  since  single-user  performance  is  achieved  in  an 
infinite-user  environment. 

In  this  chapter  we  have  analyzed  the  performance  of  nonlinear  DS/SSMA  correlation  receivers 
as  we  let  N  grow  without  bound.  Assuming  some  regularities  on  the  nonlinearity,  conditions  for 
multiple-access  signals  were  obtained  to  have  asymptotic  single-user  error  probabilities  in  an 
infinite-user  environment  This  limiting  behavior  of  Fe  implies  that  the  performance  of  multi-user 
communication  systems  in  impulsive  channels  can  ,  for  long  sequences,  be  improved  significantly, 
knowing  the  statistical  characteristic  of  the  channel  noise. 


CHAPTER  9 


SOFT-LIMITING  CORRELATION  RECEIVERS  IN  IMPULSIVE  CHANNELS: 

APPROXIMATIONS 

The  results  of  Chapter  5  indicated  that  the  improvement  in  performance  by  using  the  hard- 
limiting  receiver  in  place  of  the  linear  correlation  receivers  increases  as  the  tails  of  the  noise  distribu¬ 
tion  become  heavier  (  recall  that  the  SNR  was  fixed  ).  However,  for  less  impulsive  noise  examples 
(  e.g.,  Laplacian  noise  ),  we  observed  that  for  small  N  and  K=2  the  conventional  linear  receiver  out¬ 
performs  the  hard-limiting  correlation  receiver.  Comparing  the  two  receivers,  the  linear  correlator  is 
more  effective  against  multiple-access  interference,  whereas  the  hard-limiting  correlation  receiver  is 
more  effective  against  impulsive  channel  noise.  We  combined  these  desirable  features  and  intro¬ 
duced  soft-limiting  correlation  receivers  as  an  alternative  for  moderately  impulsive  channels.  The 
performance  of  the  soft-limiting  correlator  will  be  investigated  in  this  chapter. 

The  soft  limiter  is  the  optimum  nonlinearity  for  detecting  the  presence  of  a  class  of  signals  in 
Laplacian  noise.  In  this  case  the  memoryless  nonlinearity  is  the  clipper  function  defined  as 

{-c  x  <  — c 

x  -c  £  x  <  c  (9.1) 

c  c  <  x  , 

where  c  is  the  clipping  level.  Although  it  is  a  suboptimum  nonlinearity  in  a  multi-user  environment, 
unlike  the  hard-limiter,  the  linear  portion  of  the  soft-limiter  allows  passage  of  the  multi  pie -access  sig¬ 
nal.  Similar  to  the  hard-limiter,  the  limiting  part  of  the  soft-limiter  controls  the  effects  of  non- 
Gaussian  impulsive  noise  samples  (see  Figure  9.1). 

Computing  the  average  bit-error  probability  of  the  soft-limiting  correlation  receiver  in  a  multi¬ 
user  system  is  extremely  complicated.  Therefore,  in  this  thesis,  we  resort  to  asymptotic  performance 
analysis  and  an  upper  bound  on  the  multi-user  error  probability. 
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9.1.  Chemoff  Upper  Bound 


In  this  section,  we  bound  the  error  probability  of  the  soft-limiting  correlation  receiver  in  a  ; 

multi-user  environment  by  applying  the  Chernoff  bound.  We  begin  by  writing  the  error  probability  ' 

l 

as 


Pe(x)  =  E*,b^Pr 


Yn=0  |b^=-l,pi 


(9.2) 


1 


where  the  test  statistic  Yn  is  given  as 


N-l  *  m 

YN=  E  z/‘> 

j  =  0 


N- 1 

X  ci 


j-0 


Ti/j  +  r/l)  +  '\  yTcb^ 


(9.3) 


cl  is  the  clipper  function.  Now,  we  define  a  zero  mean  random  variable  by 


flj^Z/l)-E{Z/l)|bil>  =  -l,pI}( 


(9.4) 


where  Zj(1)  is  given  in  (9.3).  By  substituting  (9.4)  into  (9.2)  we  can  rewrite  the  error  probability  as 


Pe(l)  =  E*fb 


N-l 


z  v 

j-0 


N-l 

-  Y 

j“o 


e{z; 


(i) 


(9.5) 


Note  that  in  (9.5)  the  mean  of  Z,'n  is  the  conditional  expectation  of  Zj(I)  given  =  — I  and  is  given 

as 
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where  |i.j  =  E  { Zj  |  b^'^1  ,pi  }.  Recall  that  (9.7)  holds  for  any  positive  value  of  t;  therefore,  to 
achieve  the  best  bound  we  minimize  the  right-hand  side  of  (9.7)  with  respect  to  L  In  particular,  we 
are  looking  for  to  such  that 


to  =  arg  <  mmtm^CO  +  tjXj] 


Similar  to  the  derivation  in  Chapter  7,  differentiating  the  argument  in  (9.9)  and  setting  the  result 


equal  to  zero  yields  the  following: 


N-l  M'  (to) 

to>0- 

j-o 


(9.10) 


If  (9.10)  can  be  solved,  the  Chemoff  bound  on  the  error  probability  of  the  soft-limiting  correlation 


receiver  is  given  as 


P.CtJSE^exp 


[tom  +  m^to)]! 


(9.11) 


where  to  is  obtained  by  solving  (9.10). 

Due  to  the  complexity  of  the  expression  for  M^(t),  we  pursue  this  technique  only  for  the  Lapla- 
cian  channel  example.  The  Laplacian  density  function  is  given  by 


V’(x)  *  e 


] _ "T|x| 


(9.12) 


where  the  variance  of  r|j(1)  is  a2=  --  y  — .  The  cumulative  distribution  function  for  this  example  is 

ZEtj 
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given  as 


1  -£• 

V(a)«l--e  °  a>0. 


(9.13) 


Considering  (9.12)  and  (9.13),  the  expression  for  the  conditional  mean  of  Zj  can  be  simplified  as 


E{Zj|b<i1>=-l,pI}=-5 


V2 


1  - 1/,)  -  ^  e"  °  ^  sinh[^  (  l-lf> )  ] 


(9.14) 


Similarly,  the  moment-generating  function  of  f|j  given  in  (9.8)  can  be  simplified  to 


•JT 

1 

—  ^ 


-.cfr  1-ih 


l- 


JL 

t5a 


!+• 


M. 

t8o 


&&<52-2 


,  A  c 


(9.15) 


where  5  =  VrT'  and  c'  =  Next  (9.15)  is  substituted  into  (9.10),  which  is  then  solved  for  to- 

The  Chemoff  bound  on  the  error  probability  of  the  soft-limiting  correlation  receiver  is  obtained  via 
(9.11). 


9.2.  Asymptotic  Analysis 

In  this  section,  we  investigate  the  behavior  of  the  soft-limiting  DS/SSMA  correlation  receiver 
when  the  length  of  the  signature  sequences  used  per  data  bit  becomes  infinitely  large.  As  mentioned 
earlier  this  analysis  provides  information  about  the  limiting  behavior  of  Pe  that  could  be  useful  in 
cases  when  N  is  large  but  still  finite.  As  before,  we  apply  the  central  limit  theorem  and  obtain  a 
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Gaussian  estimate  of  the  error  probability  for  the  soft-limiting  correlation  receiver  which  is  asymptot¬ 
ically  exact 

Recall  from  the  asymptotic  analysis  in  Chapter  8  that,  conditioned  on  the  multiple-access 
parameters  pi,  the  test  statistic  is  a  sum  of  N  independent  random  variables.  Furthermore,  note  that 
the  received  data  can  be  inflated  by  VN  to  keep  the  correct  scale,  and  then  the  test  statistic  is  written 
as 


N— 1 


Yn=  E  cl 

j-0 


tlj  +  VNTj(l)  + 


(9.16) 


N0T  ... 

where  T|j’s  are  samples  with  zero  mean  and  variance  — -  -  .  In  (9.16),  are  defined  in  (8.2)  and 

(8.3).  Applying  the  central  limit  theorem  to  the  sum  (9.16),  the  asymptotic  bit-error  probability  for 
the  soft-limiting  DS/SSMA  correlation  receiver  can  be  written  as 


55  Er,$,b  |  Q 


(9.17) 


where  Ey  =  lim  E 

N— »+«• 


~  !b^  =+l,pi 


,  a 

and  (oy )  =  lim  Var 
N  — »+•• 


.  Notice  that 


(9.17)  can  be  thought  of  as  an  estimate  of  the  bit-error  probability  of  the  soft-limiting  DS/SSMA 
correlation  receiver  which  is  asymptotically  exact.  We  pursue  the  estimate  based  on  (9.17)  by  ela¬ 
borating  on  the  argument  of  the  Q  function.  For  a  given  set  of  multiple-access  parameters  the  nor¬ 
malized  mean  of  the  test  statistic  for  the  Laplacian  channel  example  is  given  as 


*V5.1  VTT  1  +  -  I- 


..  No 

where  o  —  ^p  (i>  •  Note  that  for  very  large  N,  (9.18)  can  be  closely  approximated  by 


=  [l+Ii(b,x,$)-e~N^"C[l+I1(b,T,$]  , 


(9.19) 


A  Yn 


where  Y  =  -^==.  Similarly,  the  variance  of  Y  is  given  as  follows: 


,  N0T 

Of  =  — - —  +  lim 
4  N->+ 


»* 


N-i  ““hi 

X  — 

j-0 


(1+1/°)] 


(9.20) 


where  Yj  is  the  argument  of  cl  in  (9. 16).  Also  notice  that  for  large  N,  (9.20)  can  be  approximated  by 


2  C  4  6 


(9.21) 


To  obtain  an  estimate  for  the  error  probability  of  the  soft-limiting  correlation  receiver,  (9.19)  and 
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(9.21)  are  substituted  into  (9.17)  which  yields 


[W_ 

V  No 


")P  (1) 

—  [l+I1(b,x,<D)](l~e  VH>T  ) 
N0 


1-0  + 


Apr 

V  n0t 


c)e 


(9.22) 


Note  that  the  accuracy  of  the  approximation  (9.22)  improves  with  N.  It  is  also  interesting  to  observe 
that  the  soft-limiter  analysis  in  this  chapter  reduces  to  the  analysis  for  the  linear  correlation  receiver  if 
we  let  c  =  «».  On  the  other,  hand  notice  that  when  the  clipping  level  is  set  to  zero,  the  receiver  gen¬ 
erates  zero  test  statistic,  and  the  resulting  error  probability  is  one  half.  Since  the  error  probability  of 
the  soft-limiting  correlation  receiver  is  estimated  by  a  simple  form  in  (9.22),  this  expression  can  be 
used  to  elaborate  on  the  performance  of  the  receiver.  In  particular,  (9.22)  is  compared  with  the 
expression  for  the  error  probability  of  the  linear  correlation  receiver  given  in  (6.1).  In  that  analysis,  if 

N-l 

the  distribution  of  the  sum  of  N  independent  identically  distributed  random  variables  £  T)'j  was 

j-o 

approximated  by  the  Gaussian  distribution,  the  expression  (6.1)  would  be  estimated  by 


(9.23) 


Note  that  for  large  enough  N,  (9.23)  results  in  a  good  estimate  of  the  error  probability  for  the  linear 
correlation  receiver  in  Laplacian  noise.  Now  comparing  the  approximation  for  the  performance  of 
the  soft-limiting  correlation  receiver  with  that  for  the  linear  correlation  receiver  we  notice  that  they 
differ  by  a  factor  in  the  argument  of  the  Q  function.  From  (9.22)  this  factor  is  seen  as  a  function  of 


the  clipping  level  c  and  is  given  as 


(9.24) 


Looking  closely  at  (9.24),  it  can  be  shown  that  (9.24)  is  greater  than  unity  for  all  finite  values  of  clip¬ 
ping  level  c.  Therefore,  according  to  this  application,  soft-limiting  correlation  receiver  outperforms 
the  linear  correlation  receiver  for  the  Laplacian  channel  example.  The  approximate  performance  of 
the  soft-limiting  correlation  receiver  can  be  evaluated  via  (9.22)  by  adjusting  the  signal-to-noise  ratio 
and  applying  the  technique  developed  in  Chapter  6  for  the  linear  correlation  receiver.  Particularly, 
the  adjustment  increases  the  effective  signal-to-noise  ratio  in  the  following  way: 


V3ac' 


otc  =  a 


1  -e 


N 


V2c 


/  ’ 


(9.25) 


where 


V  N„ 


9.3.  Numerical  Results 


The  numerical  results  presented  here  are  aimed  at  showing  an  improvement  in  performance  by 
using  soft-limiting  correlation  receivers  in  place  of  linear  or  hard-limiting  correlation  receivers  when 
the  channel  is  impulsive  but  not  excessively  so  (  e.g.,  Laplacian  channel ).  However,  due  to  the  com¬ 
plexity  of  analyzing  the  performance  of  nonlinear  DS/SSMA  correlation  receivers,  we  can  only  com- 


pare  the  exact  error  probabilities  of  linear  and  hard-limiting  correlation  receivers  with  the  estimates  of 
the  error  probability  of  the  soft-limiting  correlation  receiver  developed  in  the  previous  two  sections. 

The  Chemoff  upper  bound  is  obtained  by  first  substituting  (9.15)  into  (9.10)  and  then  solving 
for  to,  and  the  bound  is  then  obtained  via  (9.11).  Table  9.1  contains  results  for  the  Laplacian  channel 
example  with  signal-to-noise  ratios  8  and  4.  These  results  are  obtained  with  the  same  DS/SSMA  sys¬ 
tem  parameters  considered  in  the  examples  in  Chapter  5.  These  examples  are  carried  out  with  lengths 
of  signature  sequences  equal  to  31,  63,  127,  and  255  and  different  values  of  the  clipping  level  c. 
Comparing  Table  9.1  with  Table  7.3  the  Chemoff  upper  bound  on  the  error  probability  of  the  soft- 
limiting  correlation  receiver  results  in  lower  upper  bounds  than  that  for  the  hard-limiting  correlation 
receiver  for  all  the  values  of  N. 

Again  for  the  Laplacian  channel  example,  we  examine  the  Gaussian  approximation  introduced 
in  the  previous  section  which  estimates  the  performance  of  the  soft-limiting  correlation  receiver  via 
(9.22)  using  the  Taylor  series  expansion  technique.  These  results  which  are  asymptotically  exact  as 
N —»+«,,  are  used  in  comparing  soft-limiting  correlation  receivers  with  linear  receivers  and  are  dep¬ 
icted  in  Tables  9.2  and  9.3.  These  tables  are  generated  with  signal-to-noise  ratio  8  and  4,  respec¬ 
tively,  and  with  different  clipping  level.  The  examples  are  carried  out  with  lengths  of  signature 
sequences  equal  to  31,  63,  127,  and  255.  Here,  by  comparing  the  approximate  error  probabilities  in 
Table  9.2  with  the  exact  ones  in  Table  4.1,  considerable  improvement  in  performance  is  observed  by 
using  soft-limiting  correlation  receivers  in  place  of  linear  correlators  for  the  Laplacian  channel  exam¬ 
ple.  However,  note  that  the  average  error  probability  for  the  linear  receiver  is  an  exact  computation, 
whereas  results  for  the  soft-limiting  correlation  receivers  are  merely  approximations  of  the  average 
error  probability  and  are  only  asymptotically  exact.  Although  these  finding  are  not  decisive  they  do 
provide  a  trend  for  further  investigation  on  the  subject. 
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TABLE  9.1.  CHERNOFF  UPPER  BOUND  FOR  THE  ERROR  PROBABILITY 
OF  SOFT-LIMITING  CORRELATION  RECEIVERS  IN  THE  BINARY  PSK 
DS/SSMA  SYSTEM  ;  LAPLACIAN  CHANNEL. 


NORMALIZED  CLIPPING  LEVEL  c‘ 


SNR  =  8  dB 

31 

1.7 

4.74  x  KT4 

2.94  x  10~3 

63 

1.5 

135  xl(T* 

ZS7  x  KT4 

127 

1.5 

6.67  x  KT5 

1.24  x  1  <r* 

255 

1.4 

3.28  x  1(T5 

5.01  x  ir5 

SNR  -  A  dB 


31 

1.4 

2.67  x  10~2 

4.89  x  KT2 

63 

1.4 

1.94  x  1(T2 

2.36  xlO-2 

127 

1.4 

1.47  x  1(T2 

1.80  xl(T2 

255 

1.5 

1.21  x  1(T2 

1.39  x  10"2 

i 


t , 


TABLE  9.2.  ASYMPTOTIC  APPROXIMATION  FOR  ERROR  PROBABIL¬ 
ITY  OF  SOFT-LIMITING  CORRELATION  RECEIVERS  IN  THE  BINARY 
PSK  DS/SSMA  SYSTEM  ;  LAPLACIAN  CHANNEL,  SNR-8.0  dB. 


NORMALIZED  CLIPPING  LEVEL  c' 


1.0 

0.0521 


1.0 

0.09 


0.1 

0.148 


0.1 

0.46 


III 


4.95  x  10"7 

5.86  x  10^ 

2.60  x  10-7 

3.83  x  10-* 

3.56  x  ir7 

2.39  x  10^ 

2.53  x  1CT7 

1.89  x  10-6 

2.60  x  ir7 

6.51  x  10"7 

2.52  x  IQ-7 

6.35  x  10-7 

4.52  x  10r7 

6.88  x  10~7 

2.65  x  10 -7 

4.19  x  ir7 
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CHAPTER  10 


SUMMARY  AND  CONCLUSIONS 

In  this  thesis  we  have  considered  an  asynchronous  multi-user  communication  problem  over 
impulsive  non-Gaussian  channels.  To  have  multiple- access  capability  the  direct-sequence  spread- 
spectrum  technique  was  used.  In  the  receiving  end,  linear  and  nonlinear  correlation  receivers  were 
used  and  the  performances  of  these  receivers  were  evaluated  by  examining  their  bit-error  probabili¬ 
ties.  These  nonlinear  correlation  receivers  were  developed  by  inserting  a  nonlinear  element  in  the 
structure  of  the  conventional  linear  correlator.  This  nonlinearity  was  introduced  to  limit  the  influence 
of  the  impulsive  channel  on  the  test  statistic  used  by  the  receiver.  A  tractable  way  of  studying  corre¬ 
lation  receivers  in  non-Gaussian  impulsive  channels  was  also  introduced  by  modeling  samples  of 
noise  after  front-end  filtering.  The  main  assumption  concerning  the  additive,  zero  mean  channel 
noise  is  that  these  samples  taken  at  the  chip  rate  are  independent.  This  allows  us  to  study  impulsive 
noise  sources  by  modeling  their  post-sampling  first-order  probability  distribution  functions.  Among 
various  first-order  non-Gaussian  models  we  considered  the  e-mixture  of  two  Gaussian  distributions, 
and  the  Laplacian  distribution;  these  are  tractable,  commonly  used  empirical  models  for  impulsive 
environments.  The  multi-user  communication  system  under  study  was  examined  against  these  impul¬ 
sive  noise  channels. 

The  model  of  the  asynchronous  binary  PSK  direct-sequence  SSMA  system  considered  here 
allowed  a  number,  K,  of  users  to  share  a  channel.  The  linear  and  nonlinear  correlators  were  assumed 
to  be  matched  to  the  first  of  the  users’  signal.  The  structure  of  the  nonlinear  correlation  receiver 
included  an  integrator,  and  the  output  was  sampled  at  the  chip  rate.  These  samples  were  then  passed 
through  a  memoryless  nonlinearity.  The  decision  on  the  parity  of  a  data  bit  was  based  on  the  sum  of 
N  samples  corresponding  to  that  bit  taken  at  the  output  of  the  nonlinear  element.  Among  nonlineari¬ 
ties  we  were  primarily  interested  in  the  hard-limiter,  which  is  known  to  be  very  effective  against 
impulsive  disturbances  in  single-user  communication  problems.  In  particular,  we  compared  the 


performance  of  the  conventional  linear  correlator  with  that  of  the  hard-limiter  correlation  receiver  in 
impulsive  and  multi-user  noises.  Then  we  briefly  examined  the  soft-limiting  correlation  receiver  in  a 
similar  environment. 

To  examine  the  performance  of  linear  correlation  receivers  we  computed  the  average  bit-error 
probability  when  multi-user  and  impulsive  channel  noises  are  interfering  with  a  binary  PSK  signal. 
In  all  of  our  examples,  a  maximal-length  spreading  sequence  of  period  N  was  assigned  to  each  of  the 
K  users  to  reduce  multiple-access  interference.  The  principal  conclusion  here  was  that,  with  a  fixed 
SNR,  the  linear  correlation  receiver  does  not  perform  as  well  as  the  Gaussian  model  predicts  when 
the  non-Gaussian  noise  has  an  impulsive  nature  (  heavy-tailed  distribution  ).  In  one  case,  in  which 
the  Gaussian  assumption  was  violated  in  favor  of  an  e-mixture  assumption,  the  resulting  degradation 
in  performance  was  equivalent  to  that  caused  by  subscribing  four  additional  users  to  a  Gaussian  chan¬ 
nel  with  same  SNR. 

Motivated  by  the  poor  performance  of  the  linear  correlation  against  impulsive  noise  channels, 
we  studied  the  performance  of  the  hard-limiting  correlation  receiver  under  the  same  conditions  as 
those  for  which  the  linear  correlation  receiver  was  analyzed.  We  first  showed  that  hard-limiting 
correlation  receivers  offer  significant  improvement  in  performance  over  the  linear  correlation  receiver 
in  single-user  impulsive  channels.  Then  the  combined  effects  of  impulsive  noise  and  multi-user 
interference  on  the  performance  of  a  hard-limiting  correlation  receiver  were  compared  with  those  on 
the  performance  of  a  linear  correlation  receiver.  These  results  indicated  significant  improvement  in 
performance  by  using  hard-limiting  correlation  receivers  in  place  of  linear  correlators  for  more  impul¬ 
sive  noise  channels.  Moreover,  the  improvement  became  more  visible  as  the  length  of  the  signature 
sequences  used  by  the  channel  subscribers  increased.  We  also  noted  that,  unlike  the  linear  correlation 
receiver,  degradation  in  performance  due  to  the  interfering  user  was  no  longer  uniform  in  a  range  of 
signal-to-noise  values.  In  fact,  as  the  SNR  increases,  the  hard-limiter  apparently  is  not  as  effective  in 
separating  the  two  users  as  the  linear  correlation  is.  However,  this  deficiency  is  outweighed  by  the 
improvement  against  the  impulsive  noise  as  long  as  the  channel  noise  is  significant. 


In  the  analysis  of  hard-limiting  correlation  receivers  we  observed  that  in  examples  of 
moderately  impulsive  channel  noise  with  small  N  and  K=2  the  conventional  linear  receiver  outper¬ 
forms  the  hard-limiting  correlation  receiver.  Comparing  the  two  receivers,  the  linear  correlator  is 
more  effective  against  multiple-access  interference,  whereas  the  hard-limiting  correlation  receiver  is 
more  effective  against  impulsive  channel  noise.  We  combined  these  desirable  features  and  intro¬ 
duced  soft-limiting  correlation  receivers  as  an  alternative  for  moderately  impulsive  channels.  Due  to 
the  complexity  of  analyzing  the  performance  of  soft-limiting  correlation  receivers,  we  could  only 
compare  the  exact  error  probabilities  of  linear  and  hard-limiting  correlation  receivers  with  the  esti¬ 
mates  of  the  e-ror  probability  of  the  soft-limiting  correlation  receiver.  For  the  Laplacian  channel 
example  the  approximate  performance,  which  is  asymptotically  exact  as  N  — »+°°,  indicates  consider¬ 
able  improvement  in  performance  by  using  the  soft-limiting  correlation  receiver  in  place  of  both  the 
linear  and  hard-limiting  correlation  receivers.  Although  these  findings  are  not  conclusive,  they  do 
provide  a  trend  for  further  investigation  on  the  subject 

To  gain  further  insight  into  nonlinear  correlation  receivers  and  explore  the  fundamental  limita¬ 
tion  of  these  systems,  we  carried  out  an  asymptotic  analysis  for  linear  and  general  nonlinear 
DS/SSMA  correlation  receivers.  In  this  asymptotic  analysis  the  length  of  the  data  period,  the  energy 
per  data  bit  for  a  user,  the  noise  power,  and  in  turn  the  signal-to-noise  ratio  are  kept  constant.  How¬ 
ever,  the  length  of  the  signature  sequences  used  per  data  bit  was  allowed  to  grow  without  bound,  and 
in  turn  the  sampling  period  decreased  to  zero.  Assuming  some  regularities  on  the  nonlinearity,  condi¬ 
tions  were  obtained  under  which  a  multi-user  system  achieves  asymptotic  single-user  error  probabili¬ 
ties  in  an  infinite-user  environment  Similar  ideal  asymptotic  conditions  on  the  MA  signals  were 
obtained  when  linear  receivers  were  considered.  An  expression  was  also  given  for  the  asymptotic 
error  probability  of  linear  and  nonlinear  correlation  receivers  in  Gaussian  and  non-Gaussian  channels. 
This  limiting  behavior  of  the  probability  of  error  reinforced  the  idea  that  the  performance  of  the 
multi-user  communication  system  operating  in  impulsive  environments  can  be  improved  significantly 
by  taking  into  account  the  statistical  characteristics  of  the  channel  noise. 


APPENDIX  A 


NON-GAUSSIAN  CHANNEL  MODELS 

In  this  Appendix  we  introduce  some  classes  of  first-order  probability  distribution  functions  that 
ate  widely  used  in  modeling  channel  noise  sources.  These  classes  complement  the  e-mixture  and 
Laplacian  noise  models  discussed  in  Chapter  2. 

Among  the  physical  models  for  impulsive  noise,  some  of  the  most  general  are  those  developed 
by  Middleton  [19-20].  Canonical,  statistical-physical  models  have  been  constructed  and  experimen¬ 
tally  verified  for  a  broad  class  of  impulsive  (  mostly  man-made  )  electromagnetic  interference.  In 
addition  to  avoiding  the  limitations  of  empirical  models,  these  models  also  remain  tractable  for  most 
analyses.  Of  particular  interest  to  us  are  density  functions  in  the  Middleton  Class  A  impulsive  noise 
model  which  can  be  written  as 

V*)=  ZKn^OjC^),  (A.1) 

m=0 

where  N( 0;C^)  denotes  the  zero  mean  Gaussian  probability  density  function  with  variance  Cm-  Here 

g~A^ni  m 

Km  = - ; —  and  C,£  =  +r']  where  A,  V  and  arc  three  parameters  of  the  model.  The 

m!  A 

variance  of  a  random  variable  with  density  given  by  (A.l)  is 

02=  £  Kmc£  =  fl2A(l+n  .  (A. 2) 

m=0 

With  the  variance  of  the  random  variable  fixed,  there  are  only  two  free  parameters,  usually  taken  to 
be  A  and  T'.  The  first  free  parameter,  A,  is  called  the  impulsive  index  and  measures  the  amount  of 
temporal  overlap  among  the  waveforms  of  the  interfering  signals  ( inside  the  receiver ).  A  large  value 


of  A  means  considerable  overlap  with  a  corresponding  approach  to  Gaussianity,  while  a  small  value 
of  A  means  highly  impulsive  or  structured  interference.  The  other  variable,  Tv,  is  given  by  the  ratio 
of  the  power  in  the  Gaussian  portion  of  the  interference  to  the  power  in  the  impulsive  component 
Middleton  has  shown  that,  by  adjusting  the  parameters  A  and  T',  the  density  given  in  (A.l)  can  be 
made  to  fit  a  great  variety  of  non-Gaussian  densities  quite  well  [19-20],  and  [36-37].  The  parameters 
A  and  T*  are  physically  motivated,  unlike  the  parameters  in  the  e-mixture  model,  and  can  be 
estimated  from  physical  measurements.  (  This  problem  has  been  studied  recently  by  Zabin  [46],  ) 
Vastola  [40]  has  shown  that  in  many  cases  the  infinite  series  in  (A.l)  can  be  closely  approximated  by 
only  two  terms.  This  yields  a  way  of  choosing  parameters  of  the  mixture  direcdy  from  the  physically 
motivated  parameters  of  the  Class  A  model.  For  the  Gaussian-Gaussian  mixture,  this  relation  stands 
as 


T  =  1+- 


a  r* 


(A.3) 


and 


Ki  _  A 
Kq+K,  1+A  * 


(A.4) 


Therefore,  when  an  e-mixture  distribution  is  used  to  model  narrow  band  non-Gaussian  noise,  e  and  y2 
can  be  obtained  from  (A.3)  and  (A.4)  knowing  that  techniques  are  available  for  determining  A  and  T'. 

Among  useful  empirical  models,  the  Cauchy  density  function  is  an  example  of  a  very  heavy¬ 


tailed  distribution.  The  density  is  given  by 


127 


where  c  is  a  parameter  controlling  the  scale.  Due  to  its  extremely  heavy  tails,  the  Cauchy  density 
function  has  infinite  variance.  However,  by  carefully  choosing  the  parameter  c,  it  can  still  be  used  as 
a  model  for  comparison  with  random  variables  of  finite  variance.  This  is  done  by  equating  the  area 
under  a  truncated  Cauchy  density  function  with  that  of  a  Gaussian  density  with  the  desired  variance. 
That  is,  we  can  use  the  relationship 


bo 


l-2Q(b)=I^ 


(A.6) 


where  Q(b)  =  -^==  J  e  2  dx  is  one  minus  the  standard  error  function,  and  c2  is  the  noise  variance 


to  be  matched.  From  this  relationship  c  is  obtained  as 


c  = 


ba 


tan[y-7tQ(b)] 


(A.7) 


( In  the  sequel  we  will  choose  b=l  in  this  context.  )  The  Cauchy  density  function  with  parameter  c 
given  in  (A.7)  enables  us  to  realize  limitations  in  performance  of  the  linear  DS/SSMA  systems  in  the 
presence  of  noise  sources  exhibiting  extremely  large  impulses. 

The  last  class  of  distributions  to  be  considered  here  is  the  generalized  Gaussian  class.  In  [16], 
density  estimates  for  non-Gaussian  processes  in  the  ocean  acoustic  environment  have  been  shown  to 
agree  closely  with  generalized  Gaussian  density  for  different  exponential  decay  rates.  This  class  has 
a  symmetric  unimodal  density  obtained  by  generalizing  the  Gaussian  density  to  obtain  a  variable  rate 
of  exponential  decay  and  is  given  by  (  see  [3]  and  [22] ) 


2r(-)A(c) 

C 


where  A(c)  -  [~r~]1/2,  o2  is  the  variance,  T(.)  is  the  gamma  function,  and  c  is  a  positive  param- 

i(3/c) 

eter  controlling  the  rate  of  decay.  Note  that  for  c-2  this  density  reduces  to  the  Gaussian  density, 
whereas  for  c-1  it  becomes  the  Laplacian  ( double-exponential )  density. 


The  classes  of  density  functions  introduced  above  cover  a  wide  range  of  practical  models  for 
non-Gaussian  sources,  and  all  are  of  interest  in  practice.  We  analyze  the  performance  of  the  linear 
correlation  receiver  in  these  non-Gaussian  channels  in  Appendix  B. 
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APPENDIX  B 


PERFORMANCE  ANALYSIS  OF  THE  LINEAR  CORRELATION  RECEIVER 

IN  NON-GAUSSIAN  NOISE 


In  Chapter  3,  we  evaluated  the  single-user  average  bit-error  probability  of  the  linear  correlation 
receiver  in  e-mixture  and  Laplacian  noise  channels.  Then  in  Chapter  4,  we  added  multi-user  noise  to 
the  scenario  and  examined  the  performance  of  the  linear  receiver  in  these  muld-user  impulsive  chan¬ 
nels.  In  Appendix  A,  we  introduced  some  alternative  models  of  non-Gaussian  channels  that  were 
also  of  interest  in  practice.  In  this  Appendix,  we  examine  the  degradation  in  performance  of  the 
linear  correlation  receiver  due  to  the  channel  noise  modeled  as  Middleton  Class  A,  Cauchy,  and  gen¬ 
eralized  Gaussian  noise  in  single  and  muld-user  cases  by  computing  the  exact  average  bit-error  pro¬ 
bability. 


Recall  from  Chapter  3  that  when  K  -  1  the  average  error  probability  for  the  linear  receiver  is 
obtained  via 


1 

Pe  =  y-JT1  Ju_1(sin  u  )  <T>2(u)  du, 


(B.l) 


AT  1N 

where  02(u)  =  |E{e‘UTl0}  is  the  characteristic  function  of  the  sum  of  the  N  independent  identically 


distributed  random  variables,  rj0 ,  t\1  , . . .  with  zero  means  and  variances 


No 


2NE1P 


.  The  error  pro¬ 


bability  for  the  Middleton  Class  A  is  computed  by  first  finding  an  expression  for  the  characteristic 
function  02  ( u ).  The  characteristic  function  for  the  Middleton  Class  A  can  be  written  as 


02(11)  = 


IN 


IKme 

m«0 


(B.2) 


V 

v: 

V, 

y, 


£ 


m 


where  Km’s  and  Cm’s  are  defined  in  Appendix  A.  The  single-user  average  error  probability  is 
obtained  by  substituting  (B.2)  into  (B.l).  The  expression  for  the  <I>2 (u)  is  more  complicated  in  the 
generalized  Gaussian  example.  In  this  case  the  probability  density  function  given  in  (A.  8),  which  has 
an  exponential  form,  is  expanded  and  then  used  to  compute  the  characteristic  function.  The  final 
expression  is  as  follows 


m  F(~— — ) 

<t>2(u)  =  £  •^■u2mA(c)2m - - - 

m-0  2m!  F(l) 

c 


where  A(c)  is  defined  in  (A.8),  T(  )  is  the  gamma  function,  and  c  is  a  positive  parameter  controlling 
the  rate  of  decay.  The  error  is  then  obtained  by  substituting  (B.3)  into  (B.  1). 

For  the  Cauchy  example,  the  characteristic  function  is  as 
<D2(u)  =  e^N|u|  (B.4) 


where  c  is  obtained  from  (A.7).  Again,  substituting  (B.4)  into  (B.l)  yields  the  error  probability  for 
linear  correlators  in  Cauchy  noise,  which  can  be  written  in  a  closed  form  as 


Pe  -  — -  7t_ltan~l  (1/cN) . 


To  demonstrate  the  single-user  performance  of  linear  correlation  receivers.  Figures  B.l  and  B.2 
have  been  generated  for  these  non-Gaussian  channel  examples.  Comparing  the  non-Gaussian  channel 
examples  with  the  Gaussian  one,  these  curves  indicate  a  degradation  in  performance  for  all  impulsive 
noise  cases  over  the  entire  range  of  interest  of  signal-to-noise  ratios,  with  fairly  large  degradation  in 
some  cases.  This  is  not  surprising  since  the  linear  correlation  receiver  is  designed  to  operate  against 
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Figure  B. I.  Single-user  error  probability  for  the  linear  DS/SSMA  correlation  receiver  in  Gaussian 
Lauchy,  and  Middleton  Class  A  channels,  N-31. 
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the  Gaussian  noise  channel.  Similar  to  our  conclusion  in  Chapter  3  for  the  e-mixture  and  Laplacian 
examples,  these  figures  show  that  the  impulsive  character  of  the  channel  noise  can  drastically  degrade 
the  performance  of  conventional  linear  receivers. 

Next  we  examine  the  performance  of  the  linear  receiver  in  the  direct-sequence  multi-user 
environment  over  Middleton  Class  A,  generalized  Gaussian,  Cauchy  noise  channels.  Recall  from 
Chapter  4,  that  the  multi-user  bit-error  probability  can  be  written  as 

oo  oo 

Pe  =  rC1  Ju-1  (sin  u)  02  (u)  du  +  ft-1  Ju_1  (sin  u)  02  (u)  [l-0i(u)]  du,  (B.6) 

2  o  o 

where  02  (u)  =  [ E  {  emT,° }  ]N  as  in  (B.l)  and  0i  (u)  =  E{emI'}  with  multiple-access  interference  Ii 
given  in  (4.5).  Since  Pe  is  written  conveniently  in  terms  of  characteristic  functions,  contributions  of 
the  MA  interference  and  noise  are  distinguishable.  Note  that  the  left  part  of  the  expression  for  the 
multi-user  error  probability  has  already  been  evaluated  for  the  non-Gaussian  channels  under  study 
here.  In  particular,  the  characteristic  functions  02  (u)  are  given  in  (B.2),  (B.3),  and  (B.4)  for  Middle- 
ton  Class  A,  generalized  Gaussian,  and  Cauchy  channel  examples,  respectively.  To  evaluate  the  error 
probability  via  (B.6)  we  use  the  characteristic  function  of  the  multiple-access  interference  0i  (u) 
given  in  (4. 14).  Therefore,  the  average  bit-error  probability  of  the  linear  correlation  receiver  for  these 
non-Gaussian  channels  is  obtained  by  substituting  the  corresponding  expression  for  02  and  (4.14) 
into  (B.6).  Tables  B.l  contains  the  two-user  error  probability  for  the  Middleton  Class  A  and  Cauchy 
channel  examples  with  SNR  =  8  dB  and  N  -  31.  Corresponding  results  for  the  generalized  Gaussian 
channel  example  are  depicted  in  Table  B.2.  For  simplicity,  the  examples  are  carried  through  under 
the  assumption  that  all  users  have  the  same  power.  All  of  these  results  are  obtained  with  fixed 
signal-to-noise  ratio,  thereby  showing  the  effects  of  the  noise  density  shape  on  performance. 

The  principal  conclusion  here  is  that  the  linear  receiver  does  not  perform  as  well  as  the  Gaus¬ 
sian  model  predicts  when  the  non-Gaussian  noise  has  an  impulsive  nature  (  heavy-tailed  ).  For 


TABLE  B.l.  AVERAGE  ERROR  PROBABILITY  OF  LINEAR  CORRELA¬ 
TION  RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ;  MIDDLE- 
TON  CLASS  A,  CAUCHY  AND  GAUSSIAN  CHANNELS,  SNR-8.0  dB, 
N-31 


DISTRIBUTIONS 

K 

i  waammmm 

MIDDLETON  CLASS  A 

1 

4.59  x  10-4 

A  -  0.0001,  r-  50.0 

2 

6.72  xl( r* 

1 

4.80  x  1(T* 

A  -  0.35,  r  -  0.0005 

2 

7.72  x  10^ 

I 

1.04  x  10-3 

A-o.i,r-o.i 

2 

1.41  x  10-3 

6 

4.66  x  ir3 

1 

1.91  x  10~* 

GAUSSIAN 

2 

4.16  x  10^ 

6 

3.48  x  10*3 

1 

2.25  x  10"1 

CAUCHY 

2 

2.26  x  lOr1 

6 

3.30  xir1 
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TABLE  B.2.  AVERAGE  ERROR  PROBABILITY  OF  LINEAR  CORRELA¬ 
TION  RECEIVERS  IN  THE  BINARY  PSK  DS/SSMA  SYSTEM  ;  GENERAL¬ 
IZED  GAUSSIAN  CHANNELS,  SNR-8.0  dB,  N-3 1. 


DECAY  RATE  c 

K 

P. 

1 

2.94  x  10-4 

1.0 

2 

5.45  x  KT4 

1 

1.91  x  10*4 

2.0 

2 

4.16  x  KT4 

1 

1.72  xltr4 

3.0 

2 

3.91  x  l(r* 

1 

1.64  xlCT4 

4.0 


2 


3.82  x  10T* 
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